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ABSTRACT 


THE  NATURE  OF  BXKLING  IN  THIN  SPHERICAL  SHELLS 
by  Lynn  Seaman 


The  report  deals  with  a  number  of  questions  concerning  the  buckling 
of  thin  elastic  spherical  shell  segments  under  uniform  external  pressure. 
Both  experimental  and  analytical  studios  were  performed  to  determine  the 
nature  of  the  buckling  phenomenon. 

The  experimental  work  was  performed  on  plastic  shells  vdiich  buck¬ 
led  elastically  so  that  tests  coxild  be  repeated.  The  size  and  position 
of  the  buckle  or  dimple,  effect  of  creep  on  buckling  loads,  importance 
of  the  precision  of  clamping  conditions,  and  the  reproducibility  of 
buckling  loads  under  unchanged  conditions  were  determined.  A  general 
discussion  was  given  of  the  shape  of  the  load-deflection  curve,  the 
appearance  of  the  buckle  through  the  post-buckling  range  and  the 
dependence  of  buckling  pressures,  unbuckling  pressures  and  minimum 
pressxires  on  the  shape  parameter  ^  .  No  correlation  was  possible 
between  the  observed  imperfections  and  the  buckling  loads. 

In  the  analysis,  the  nonlinear  equation  for  total  potential 
energy  was  derived  for  deep  shell  investigations.  A  formula  for 
deflected  shape,  which  was  similar  to  that  seen  in  the  tests,  was 
used  to  determine  equilibrium  positions  for  a  number  of  shell  con¬ 
figurations.  Besides  the  usual  shape  parameter  >^  ,  a  depth  measure 
was  also  necessary  since  the  shells  were  not  shallow.  The  equililadum 
positions  were  plotted  into  load-center  deflection,  load-average  de¬ 
flection,  and  energy-deflection  curves.  Computed  buckling  loads  were 
much  too  high,  probably  due  to  neglecting  the  mode  shapes  which  are 
small  surface  undulations  which  occur  in  the  shell  before  buckling. 

The  equation  developed  is  adequate  in  the  range  near  buckling  but  needs 
improvement  for  use  in  the  largo  deflection  range. 

The  experimental  buckling  loads  were  in  a  fairly  narrow  band  midst 
those  of  other  investigators.  The  repeatability  of  the  tests  added 
greatly  to  the  value  of  the  results.  Analytical  buckling  loads  wore 
higher  than  previous  resvilts,  and  far  above  the  test  results. 
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CHAPTER  1  INTRODUCTION 


With  the  increased  use  of  shell  type  structures  in  architecture, 
space  exploration,  and  other  projects,  much  interest  has  been  evinced 
In  the  determination  of  the  collapse  or  buckling  loads  of  these 
structures.  This  report  treats  the  buckling  of  thin  elastic  spherical 
shell  segments  under  the  action  of  uniform  external  pressure. 

The  object  was  to  define  the  nature  of  the  buckling  phenomenon. 

This  Included  the  buckling  load,  shape  of  the  load -deflection  curve, 

the  effect  of  imperfections  and  testing  conditions,  and  the  size  and 

shape  of  the  dimple  which  forms.  Besides  the  buckling  load  ), 

other  pressures  of  interest  were  the  load  to  which  the  shell  Jumped 

at  buckling  (qj),  the  load  to  which  It  Jumped  at  unbuckling 

and  the  minimum  post  buckling  pressure  (q  .  ).  A  further  purpose 

nin 

was  to  find  out  If  there  was  a  difference  between  buckling  loads 
obtained  from  tests  with  displacement  control  and  those  obtained  from 
tests  with  pressure  control. 

In  the  analytical  £«tudy  the  object  was  to  derive  equillbriim 
eqpations  for  the  symmetrical  deformations  of  deep  shells  and  to  use 
these  equations  for  buckling.  The  attempt  was  to  find  a  deflected 
shape  andbuckllng  loads  and  to  determine  whether  the  equations  were 
adequate.  With  a  series  of  equilibrium  positions,  load-deflection 
curves  could  be  constructed  and  compared  with  test  results. 

The  relevant  literature  was  reviewed  in  search  of  clues  to  the 
most  successful  approach  to  the  problem.  G.  A.  Thurston,  B.  Budiansky, 


and  others  have  investigated  analytically  the  shape  of  the  load-center 
deflection  and/or  load-average  deflection  curves  and  buckling  loads. 

B.  Budiansky  and  W.  L.  Chen  have  initiated  the  discovery  of  the 
importance  of  imperfections  on  buckling  loads.  A.  Kaplan  and  Y.  C.  Fung 
have  provided  experimental  knowledge  about  deflected  shape  as  well  as 
buckling  loads.  T.  von  Karman  and  H.  S.  Tsien  have  suggested  energy 
criteria  for  buckling  which  contrast  with  a  definition  of  buckling 
pressure  as  x,hc  first  maximum  point  on  the  load-displacement  curve. 

The  experiments  were  conducted  on  thin  plastic  shell  segments  with 
depth  to  span  ratios  from  l/SOO  to  l/ll.  For  the  shells  the 
shape  parameter,  ranged  from  3.5  to  25.  A  function  of  both  the  thick 
ness  to  radius  ratio  and  shell  depth,  this  parameter  has  been  found  to 
govern  all  shallow  shell  buckling  phenomena. 

'/^  S/ff 

where  a  is  the  shell  radius,  h  is  thickness,  and  CX  is  one  half  the 
opening  angle  of  the  shell. 

Jiost  of  the  tests  were  conducted  by  displacement  control  using 

water  pressure  above  and  below  the  shell.  These  tests  provided  a 

complete  load-deflection  curve  from  no  load  into  the  post-buckling  range 

and  back  to  zero  load.  A  few  tests  were  run  with  an  approximation  to 

pressure  control.  In  all  cases  the  shells  behaved  elastically  and  the 

tests  could  be  repeated.  Besides  those  topics  mentioned  above,  the 

study  included  the  effect  of  buckle  position  and  creep  on  the  buckling 

load  and  the  variation  of  the  pressures  q  ,  Qo .  <3  ,  and  a  with  A. 

cr  ^  un  ‘mm 
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The  results  were  used  to  check  the  possible  validity  of  the  energy 
criterion.  In  all  this  study  the  reproducibility  of  the  tests  was 
very  helpful,  permitting  even  an  estimate  of  the  Importance  of  effects 
which  have  caused  so  much  scatter  in  test  results. 

Analytically,  equilibrium  positions  were  determined  by  the  principle 
of  minimum  potential  energy.  The  energy  equation  wai  derived  with  the 
aid  of  an  order  of  magnitude  analysis  and  applied  to  finding  buckling 
loads  for  a  number  of  configurations  like  the  shells  tested. 

The  experimental  work  is  presented  in  Chapter  4  and  the  analytical 
study  in  Chapter  5.  The  results  from  these  two  portions  are  brought 
together  in  the  Summary,  Chapter  6,  and  compared. 
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CHAPTER  2  NOTATION 


shell  radius 


a 

b 


E 

El 


angle  defining  edge  of  shell,  see  Fig,  1 
variable  defining  buckle  extent 

ratio  of  extensiona]  to  bending 

stiffness 

linear  and  nonlinear  components  of  deflection  used 
in  an  example 
modulus  of  elasticity 
modulus  at  one  minute  after  loading 

EaoiEii >EioiEoi <^21 >^12 1^02 >^03  >^04 

coefficients  of  energy  terms  defined  by  Eq.  (5-28) 

£  strain 

strain  at  buckling 

strain  in  6,  Z  directions  or  X,  Y,  Z  directions 
shear  strain  between  0  and  Z  directions  or  X  and  Z 
directions 

shear  strain  between  Y  and  Z,  X  and  Y  directions 
strains  in  the  -liddle  surface  in  ji,  e  directions 

®11>  ®23>  ®i3»  ®33 

linear  strains  corresponding  to  )  ^33 

g  creep  coefficient 

h  thickness  of  shell 


'cr 


43 


f  B 
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central  height  of  shell 


H 

0 

(D 

kj ,  kj 
K(t,g) 

K(T.g) 


L 

m 


nondimensional  curvatures  of  the  middle  surface  in 
the  e  directions  (curvature  times  radius) 
position  coordinate,  see  Fig.  1 

variation  of  tangential  displacement  through  thickness 
stiffness  coefficients  used  in  an  exsunple 
a  function  of  creep  and  the  time  since  commencement 
of  loading 

a  function  of  creep  and  duration  of  loading;  it 

defines  the  effect  of  creep  on  the  buckling  load 
S/rtdi  =  S/'n  >  shell  shape 

parameter 

dimensionless  volume  change,  Eq,  (5-26) 
quantity  assumed  to  be  greater  than  10  used  to  define 
orders  of  magnitude 
a  constraint 

Poisson's  ratio,  taken  as  0.3  In  the  analysis 
rotation  in  ^-Z  plane 

nondimensional  pressure 
critical  or  buckling  pressure 

nondimensional  pressure  to  which  shell  jumps  at 
buckling  and  at  unbuckling,  respect.. /ely 
normalized  pressure,  see  Eq.  (4-4) 
total  potential  energy 
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fr 


*^min’  '^max 


q 


cr 


q 


m 


qi 


Qa 

Q 


S 


(T 

^1  >  ^  >  ^3 
<^r 

t 


t 

O 

r 

T 

u 

c 

u 


dimensionless  form  of  potential  energy 
pressure 

linear  buckling  pressure 

pressures  of  minimum  and  maximum  points  on  the  load- 
defiection  curve 
buckling  pressure 

definition  of  buckling  pressure  given  by  K.  0.  Friedrichs 

upper  buckling  pressure  as  defined  by  H.  S.  Tsien 

pressure  to  which  shell  jumps  at  buckling 

normalized  buckling  pressure,  see  Eq.  (4-4) 

dimensionless  pressure,  see  Eq.  (5-26) 

radius  of  shell  segment  in  plan 

shell  surface 

stress 

stress  in  ©,  Z  directions 
critical  stress 
normalized  critical  stress 

time  from  beginning  of  loading 

--i 
e  ® 

variable  of  integration  used  in  Eq.  (4-3) 
time  from  beginning  of  loading  to  buckling 
deflection  tangential  to  middle  surface 
nondimensional  deflection  of  middle  surface  tangen¬ 
tial  to  the  middle  surface  (def lection/radius ) 
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u 


strain  energy 


Ui,  Uj,  U3, 

i 

V 

V 

w 

c 

w 


w 

o 


W 


X,  Y,  Z 
X,  Y,  Z 
Z 


extensional  and  bending  conponents  of  strain  energy 

a  dimensionless  strain  energy  terms,  see  Eq.  (5-26) 
b 

position  coordinate,  see  Fig.  1 
displacement  component  in  Y  direction 
volume 

deflection  normal  to  the  middle  surface 
nondimensional  deflection  of  the  middle  surface  normal 
to  the  middle  surface  (def lection/radius ) ,  or 
displacement  in  the  Z  direction 
nondimensional  linear  deflection 
nondimensional  depth  of  buckle  at  apex  of  shell 
(depth/radius ) 
work  of  the  external  forces 

variation  of  normal  displacement  through  the  shell 
thickness 

rectangular  coordinates 

surface  loads  in  the  coordinate  directions 
position  coordinate,  see  Fig.  1 
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CHAPTER  3  HISTORY  OF  THE  PROBLEM 


The  first  attempts  to  solve  the  buckling  problem  dealt  with  the 

linearized  equations  and  assumed  that  buckling  occurred  when  two 

adjacent  equilibrium  positions  could  be  found  for  the  same  load.  The 

solution  of  the  case  with  axisyi.jnetric  deformations  was  given  by 

R.  Zoelly^^^  in  1915.  This  derivation  can  be  seen  in  Timoshenko:  The 

(2} 

Theory  of  Elastic  Stability^  .  The  critical  pressure  is  given  as 


A.  Van  der  Neut^  later  showed  that  unsymmetric  buckling  modes  did 

not  yield  a  lower  critical  value.  Hence  the  linear  buckling  load  for 

a  complete  spherical  shell  was  established  and  served  to  define  the 

important  parameters  in  the  problem.  However,  experimental  buckling 

loads  were  generally  ^  to  i  of  this  linear  value  so  other  criteria  and 

methods  of  analyses  were  tried  on  the  problem. 

An  important  idea  was  the  suggestion  that  the  shell  could  jump 

fiom  one  position  on  the  load-displacement  curve  to  a  position  with 

lower  or  equal  total  potential  energy  if  one  were  available.  (A 

possible  curve  of  pressure  versus  deflection  is  shown  in  Fig.  2.) 

(4) 

T.  von  Karman  and  H.  S.  Tsiea  announced  the  notion  of  an  energy 

(5  ) 

criterion  b\  defining  q  =  q  .  .  For  the  critical  load,  K.  0.  Friedrichs 
'  cr  h'ln 

suggested  q^,  the  pressure  at  a  point  on  the  initial  part  of  the  curve 
which  had  the  same  energy  level  as  a  point  on  the  post  buckling  curve. 
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T.  von  Karman  and  H.  S.  Tsien,  K.  M.  Mushtari 


,  and  V .  I . 


(6) 

(7  ) 

Feodosev'  were  investigators  who  used  the  energy  or  minimum  load 

criterion  that  q  =  q  ^  ,  Friedrichs  used  his  criterion  q  =  q  . 

cr  ^min  cr  ’m 

All  five  used  an  energy  method  of  analysis  and  all  except  Friedrichs 

analyzed  a  small  segment  of  the  dome  as  an  isolated  piece.  Von  Karman 

and  Tsien,  the  initial  investigators,  assumed  that  £^2  =  ^ong 

other  things,  but  found  a  of  .316  q^  ,  a  goou  estimate  of  buckling 

pressures.  This  investigation  also  successfully  predicted  the  size 

of  the  buckle.  Friedrichs  discarded  many  of  the  assumptions  of  von 

Karman  and  Tsien,  refined  the  equations  of  equilibrium,  and  pointed 

out  certain  handy  mathematical  techniques  which  can  be  used  on  the 

problem.  His  equations  indicated  that  q  .  could  be  negative,  thus 

min 

showing  the  falsity  of  either  his  solution  or  the  energy  criterion. 

It  is  suggested  here  that  his  neglect  of  the  boundary  conditions  caused 
occurrence  of  a  negative  load. 

The  first  investigators  were  bothered  by  the  edge  conditions 

on  the  small  segment  they  analyzed.  Feodosev  overcame  this  problem 

by  specifying  a  deflected  shape  which  cculd  match  the  slope,  shear, 

and  moment  induced  in  the  rest  of  the  shell.  He  found  a  negative 

value  for  q_.  . 

^in 

M.'ishtari  ,  using  different  equations  for  deflected  shape,  suid 
properly  applying  boundary  conditions,  found  a  small  positive  value 
for  (28^  and  16^  of  q,j,)  in  his  two  analyses.  He  also  critically 

examined  the  analyses  of  von  Karman  and  Tsien,  Friedrichs,  auid  Feodosev 
and  pointed  out  errors  in  their  calculations. 
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(8) 

In  a  separate  paper  following  Friedrich's  publication,  H.  S.  Tsien 
augmented  the  theories  of  von  Karman  and  Friedrichs  by  proposing  that 

the  shell  should  follow  an  energy  defection  curve,  rather  than  a  load- 

deflection  curve.  He  showed  tnat  a  negative  value  of  would  be  obtained 

from  the  von  Karman  criterion  and  that  an  additional  restriction  had  to 

be  made.  Beside  tie  condition  that  the  totai  potential  energy  must  be 

the  s8ime  before  and  after  the  jump,  he  required  that  either  the  load 

must  be  the  same  (constant  pressure  case)  or  the  volume  enclosed  under 

the  shell  must  not  change  (constant  volume  case).  In  his  analyses  Tsien 

obtained  different  critical  ’oads  in  the  two  cases,  q  for  the  constant 

cr 

pressure  case  being  about  half  that  for  constant  volume.  A  number  of 
tests  performed  under  his  direction  verified  his  results  for  constant 
volume  and  one  test  gave  some  verification  for  his  results  for  constant 
pressure . 

The  energy  criteria  have  generally  been  applied  only  to  deep  shells 
or  to  shells  with  o<  at  least  4^  or  greater.  The  reason  for  this  is  that 
the  analyses  have  been  performed  on  the  region  of  the  shell  where  the 
dimple  should  form,  neglecting  the  rest  of  the  shel  ,  and  the  total 
angle  of  the  dimple  is  7°  or  8° 

These  criteria  predict  some  possible  load  deflection  path  which 
IS  lower  than  the  route  over  the  luaximum  point  of  the  curve.  The 
peak  of  the  curve  and  thus  the  buckling  load  may  be  reduced  by 
imperfections  or  clamping  conditions.  Since  the  energy  criteria  can 
be  used  to  find  a  minimum  possible  buckling  load,  the  criteria  should 
be  useful  in  predicting  the  effect  of  imperfect  conditions. 


16 


The  energy  criteria  and  the  analyses  based  on  them  lead  to  a 
number  of  effects,  some  of  which  can  be  verified  experimentally. 

The  criteria  require  the  presence  of  adequate  external  energy  to 
produce  a  jump  and  thus  it  seems  that  buckling  could  occur  at  a  range 
of  pressures  and  not  be  very  reproducible.  Presumably  the  critical 
loads  given  in  the  analyses  were  minimum  values  at  which  '^uckllng 
could  occurand  therefore  test  results  should  be  equal  to  or  above 
these  theoretical  values.  Th^  difference  in  buckling  loads  predicted 
for  constant  volume  and  constant  pressure  tests  is  probably  the  most 
significant  effect  and  the  most  readily  verified.  It  is  significant 
because  in  practical  situations  of  buckling  under  wind  or  snow  loading, 
the  failure  is  with  nearly  constant  pressure  conditions.  Yet  the  tests 
of  Kaplan  and  Fung^^^,  Homewood,  Brine  and  Johnson and  most  of  the 
tests  of  Tsien  were  made  with  volume  control.  If  it  is  a  fact  that 
buckling  pressures  with  constant  pressure  tests  are  much  lower  than 
for  constant  volume  conditions,  this  is  important  to  know. 

A  number  of  investigators  have  attempted  to  solve  the  differential 
equations  of  equilibrium  for  the  shell  segment  and  have  usually  considered 
that  the  buckling  load  was  the  maximum  point  on  the  load-deflection  curve. 
Cylindrical  coordinates  have  been  employed  and  the  analysis  has  been 
restricted  to  fairly  shallow  shells  (rise  to  span  less  than  l/S).  The 
technique  used  has  been  to  expand  the  dependent  variables  in  terms  of 
a  power  series  of  the  center  deflection  to  thickness  ratio.  Then  the 
coefficients  of  the  series  are  required  to  satisfy  the  governing 
differential  equations.  A.  Kaplan  and  Y.  C.  Fung^^^,  R.  R.  Archer^^^\ 
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^  X  2  3l  ' 

and  H.  J.  Weinxtschke^  '  have  used  variations  on  this  method  and  found 

it  gives  results  near  test  values  only  for  very  shallow  shells.  It  can 

be  pointed  out  from  the  analytical  results  of  Weinitschke  and  Archer 

that,  with  enough  terms  in  the  series,  an  accurate  value  for  critical 

load  can  be  found.  (That  is,  accurate  theoretically,  but  not  necessarily 

near  experimental  results)  Therefore,  it  is  apparent  that  the  deflected 

shape  is  a  series  in  which  an  extremely  large  number  of  terms  is  importan 

The  energy  or  variational  approach  has  been  used  with  the  deflection 

expanded  in  a  power  seiies.  With  this  method  the  investigators  have 

(13  ) 

found  results  only  up  to  X  =  7.  G.  P.  R.  von  Willich  derived  the 

necessary  equations  and  found  results  using  only  one  parameter. 

(14  ) 

W.  L.  Chen'  used  the  same  technique  and  two  parameters.  A.  B.  Caseman 
extended  the  investigation  to  three  parameters. 

Chen  also  investigated  the  effect  due  to  various  types  of  symmetric 

imperfections  such  as  grooves  at  midheight  in  the  shell.  His  results 

show  that  the  effect  of  an  imperfection  may  be  either  positive  or 

negative  depending  on  the  magnitude  of  the  imprefection  and  the  depth 

of  the  shell.  Hence  a  correlation  between  imperfections  and  buckling 

loads  is  rather  complex. 

In  dealing  with  shallow  shells,  many  investigators  have  replaced 
the  spherical  segment  by  a  paraboloid  with  the  same  height  and  plan 
dimensions.  The  approximation  has  the  effect  of  introducing  a  symmetric 
imperfection  with  a  magnitude  of  about  .002(h*/a)X*  at  midheight  on 
the  shell.  It  can  be  seen  from  the  formula  that  the  error  increases 
rapidly  with  7^ . 


t. 


(15) 
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H.  B.  Keller  and  E.  L.  Reiss^^^^^  have  used  a  finite  difference 
approach  fairly  successfully  in  the  shallow  shell  range  to  define 
They  also  found  qj  and  q^,  the  intermediate  value  suggested  by 
Friedrichs . 

Four  recent  attempts  have  been  eminently  successful  in  determining 

buckling  loads  for  a  large  range  of  A.  These  are  of  particular  impor- 

(12b  ) 

tance  since  they  essentially  agree  on  critical  loads.  Weinitschke 

augfmented  his  earlier  power  series  approach  by  expanding  the  series  from 

both  the  center  and  edge  and  matching  the  terms  at  an  interior  point. 

This  decreased  the  convergence  problem  a  great  deal  because  both  series 

were  shorter  than  a  single  series  expanded  about  the  apex,  and  the 

(17  ) 

accuracy  was  correspondingly  higher.  G.  A.  Thurston'  used  a  numerical 

integration  technique  to  find  solutions  at  various  points  up  to  9. 

(18  ) 

B.  Budiansky'  '  transformed  the  governing  differential  equations  into 
two  simultaneous  integral  equations  and  found  good  results  up  to /Is  13. 
With  the  energy  method  Casem^'n  found  results  up  to  7.  Since  these 
four  methods  are  essentially  different  from  each  other,  yet  their  results 
agree  very  well,  it  is  reasonable  to  conclude  that  the  buckling  loads 
for  the  perfect  shell  under  symmetric  deformation  have  been  determined 
in  the  shallow  range.  The  buckling  loads  found  in  these  analyses  are 
plotted  in  Fig.  11. 

It  may  be  mentioned  that  all  the  deflected  shapes  used  in  the 
preceding  analyses  have  been  based  purely  on  mathematical  considerations, 
not  on  the  appearance  of  the  shell  surface  during  loading. 
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Budiansky  also  considered  the  possibility  of  cesslrai  imperfuc lions 
and  found  a  significant  decrease  in  load  for  large  imperfections.  The 
central  imperfection  is  deceptive  since  it  implies  a  change  in  radius 
as  well  as  an  axisymmetric  groove  at  midheight.  If  the  central  height 
is  assumed  to  determine  radius,  then  the  imperfections  of  Chen  and 
Budiansky  are  found  to  be  quite  similar.  In  fact,  b,  the  imperfection 
magnitude  of  Chen  is  3/l6(^H),  where  is  the  error  in  central  height 
in  Budiansky 's  notation.  The  effect  of  simply  an  error  in  height 
(and  therefore  in  radius)  can  be  found  from  the  linear  buckling  load, 
rewritten  as  follows 


q 


T 


eh" 


where  r  is  the  radius  of  the  shell  in  plan  and  H  is  the  central  height 
of  the  shell  segment.  A  small  change,  AH,  in  H  gives  a  corresponding 
change  in  of 


Therefore 


Eh"" 


Mr  -  _  :?(£H)  .  2£  (^-2) 

?r  '  H 

and  is  the  nondimensional  error  in  central  height  according  to 
Budiansky.  He  obtained  reductions  of  10^  to  30^  in  buckling  load  for 
£=0.05.  A  decrease  of  10^  caused  merely  by  the  change  in  radius 
could  be  found  from  Eq .  (3-2).  In  addition,  q  is  a  function  of 
as  well  as  of  the  variables  in  q,^  and  ^  in  turn  is  a  function  of 
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versus 


yj~^.  Therefore  a  change  in  H  should  shift  the  curve  of 

Aff  U 

^  (see  Fig,  11)  down  by  q  ^  and  to  the  right  by  i  —  A  ■  This 

predictable  effect  seems  to  be  verified  by  the  plots  of  buckling 

load  given  by  Budiansky  {and  not  reproduced  herein). 

Budiansky  and  Weinitschke  have  suggested  that  probably  the 

consideration  of  unsymmetrical  deformations  will  help  to  provide  critical 

pressures  closer  to  experimental  values.  Weinitschke  stated  in  reference 

12b  that  the  circumferential  stresses  which  he  obtained  were  sometimes 

as  high  as  twice  the  meridional  stresses,  thus  suggesting  the  possibility 

(19 ) 

of  nonaxisymmetric  buckling.  A  simple  linear  analysis  by  Rabotnov 
showed  that  buckling  should  occur  when  either  of  the  two  stresses  reached 
a  critical  value. 

Extensive  sets  of  tests  have  been  performed  by  Xaplan  and  Fung; 

(20) 

Tsien;  Homewood,  Brine,  and  Johnson;  and  K.  Kloppel  and  0,  Jungbluth ' 

The  u'lckling  loads  obtained  by  these  investigators  are  shown  in  Fig. 

12.  Kaplan  and  Fung  performed  an  essential  service  by  carefully  testing 
a  number  of  magnesium  shells  with  A  values  between  4  and  10.  They  not 
only  found  buckling  loads,  but  also  measured  initial  imperfections  ana 
deflected  shapes  as  the  shell  approached  buckling.  The  measurements 
were  plotted  to  show  that  during  loading  (before  buckling)  the  shell 
could  form  either  a  central  dip,  or  a  jyinmetric  furrow  about  the  shell 
between  center  and  edge,  or  both  dip  and  furrow,  or  5 jme  unsymraetric 
shape.  Furthermore,  the  particular  form  depended  on  the  value  of  Af 
the  more  complex  shapes  occurring  with  the  larger  A- 
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These  different  deflection  patterns  have  usually  been  referred  to 
as  modes  of  deflection  and  appear  quite  similar  to  vibration  modes. 

Reiss,  Greenberg,  and  Keller Weinitschke,  and  Budiansky  have  also 
mentioned  the  changes  of  mode  shapes.  As  A  is  varied  the  mode  seems 
to  change  at  about  3.2,  5.3,  8.8,  16,  25,  .  .  .  (These  are  equal 
logarithmic  steps).  The  first  value  of  A  is  about  the  minimum  at  which 
buckling  occurs.  Between  3.2  and  5.3  a  central  dip  forms  during  loading. 
Between  5.3  and  8.8  the  maximum  deflection  may  occur  in  a  symmetric 
furrow  about  the  shell.  Beyond  this  point  there  is  confusion  as  to 
whether  the  mode  is  symmetric  or  unsymmetric,  and  the  changeover 
points  are  not  well  defined. 

The  tests  of  Tsien  were  for  shells  in  a  deeper  range,  from^=  15 

to  25  and  the  buckling  loads  were  equal  to  or  slightly  a±>ove  those  of 

Kaplan  and  Fung . 

Kloppel  and  Jungbluth  tested  a  large  number  of  ribbed  domes  to 
find  the  effect  of  stiffeners  on  buckling  strength.  Their  shells 
were  made  by  welding  a  series  of  meridional  and  circumferential 
ribs  to  a  thin  shell  dome  so  that  the  structure  became  a  composite  of 

thin  shell  and  frame  types.  There  is  no  theoretically  predicted 

buckling  load  for  such  a  structure  so  it  was  assumed  to  be  a  uniform 
shell  with  some  effective  thickness. 

An  expression  which  can  be  used  for  this  effective  thickness  is 

Sg-Sg  (3-3) 
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where  S  is  an  average  extensional  stiffness  (Kh  for  a  uniform  shell) 

£ 

and  Sg  is  the  average  bending  stiffness  (-j^ch), 

Kloppel  and  Jungbluth  tried  many  combinations  of  rib  sizes  and 
distributions  yet  their  buckling  loads  were  rather  close  together  as 
seen  in  Figure  12.  Their  procedure  for  finding  an  effective  thickness 
(which  is  similar  to  Eq .  3-3)  was  somewhat  validated  by  this  lack  of 
dispersion.  Tests  were  performed  under  air  pressure  so  that  when  failure 
occurred  it  was  sudden  and  complete.  Many  of  the  specimens  exhibited 
local  failure  before  the  snap-through  tock  place.  It  may  be  postulated 
that  local  failure  occurred  in  all  cases,  but  it  was  only  observed  in 
those  where  buckling  did  not  occur  immediately  thereafter. 

The  buckling  loads  of  Kloppel  and  Jungbluth  are  considerably  lower 
than  those  of  Tsisn  in  the  same  range  and  this  fact  may  be  explained 
partially  by  the  occurrence  of  local  buckling  and  partially  by  the 
assumed  equivalence  of  frame  and  uniform  shells.  Since  many  dome 
structures  are  actually  composite,  the  work  of  these  two  investigators 
has  a  great  deal  of  practical  significance. 

Quite  recently  R.  H.  Homewood,  A.  C.  Brine,  and  A.  E.  Johnson 
made  a  series  of  tests  with  ^  from  6  to  20,  thus  connecting  the 
experimental  ranges  of  Kaplan  and  Fung  and  Tsien.  Their  results  agreed 
well  with  Kaplan  and  Fung  but  were  considerably  below  the  buckling 
loads  of  Tsien,  Two  of  their  shells  unbuckled  when  the  pressure  was 
removed  showing  that  the  material  had  not  been  stressed  beyond  the 
elastic  range.  These  two  were  subsequently  re-buckled  and  the  second 
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buckling  loads  ware  nearly  equal  to  the  critical  loads  obtained  the 
first  time.  Hence  some  reproducibility  of  buckling  loads  could  be 
deduced  from  these  tests. 

The  combined  test  results  show  a  fairly  large  scatter,  particularly 
between  buckling  'oads  of  different  investigators.  The  dispersion  between 
test  results  of  separate  experimenters  may  be  attributed  to  differences 
tn  shell  material,  test  apparatus  and  procedures,  clamping  conditions, 
and  the  range  of  A  used. 

The  analyses  and  tests  performed  tnus  far  have  been  adequate  to 
solve  a  few  problems  of  buckling  and  to  raise  a  myriad  of  questions. 

From  the  tests  It  can  be  predicted  that  q^^  is  between  15^  and  80^ 
of  the  linear  buckling  pressure.  The  shell  surface  does  not  remain 
uniform  during  loading  but  forms  waves  and  the  number  of  waves  is 
partially  predictable.  A  curve  of  buckling  pressure  versus  A  has 
been  found  analytically  for  symmetric  deformations  of  shallow  shells, 
but  this  pressure  does  not  correspond  well  with  test  results.  The 
convergence  of  iteration,  perturbation,  and  finite  difference  techniques 
often  is  not  adequate  near  the  buckling  point.  This  suggests  that  the 
correct  deflection  shape  ’s  a  power  series  with  a  large  number  of 
important  terms.  Imperfections  have  been  considered  both  analytically 
and  experimentally  and  so  far  the  only  indication  is  that  the  effect 
ot  imperfections  la  rather  complex. 


24 


The  scatter  in  test  results  and  the  disagreement  between 


theoretical  and  experimental  buckling  loads  prompted  the  author  to 
consider  a  number  of  questions.  The  investigations  of  this  report  were 
designed  to  provide  some  answers  for  these  problems.  Some  of  the 
pertinent  questions  were  the  following: 

1.  Are  imperfections  important?  What  is  the  effect  of  size, 
shape  and  position  of  imperfections  on  buckling? 

2.  Are  clamping  conditions  important? 

3.  Are  buckling  deformations  symmetric? 

4.  Are  buckling  loads  reproducible? 

5.  Is  an  energy  criterion  valid  so  that  the  critical  loads 

depend  on  the  ambient  energy  level? 

6.  What  is  the  position  of  the  dimple  which  forms  at  buckling 
and  does  the  position  affect  the  critical  load? 

7.  Does  the  type  of  loading  —  constant  volume  or  constant 
pressure  —  affect  the  buckling  pressure? 

8.  Is  buckling  precipitated  by  yielding  or  local  failure  in 
cases  where  q^^  is  below  the  analytically  predicted  pressure? 
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CHAPTER  4  EXPERIMENTAL  STUDY 


A  Introduction 

The  purpose  of  this  experimental  project  was  to  investigate 
the  nature  of  buckling  of  thin,  clamped,  spherical  shells  under  uniform 
pressure.  Under  this  general  topic  the  following  specific  items  were 
considered : 

a.  Shape  of  the  load-displacement  curve  from  no  load  to  post- 
buckling  and  the  rebound  or  unbuckling  curve. 

b.  Size,  shape,  and  position  of  the  buckle  which  formed. 

c.  Reproducibility  of  buckling  loads. 

d.  Effect  of  shell  imperfections  or  clamping  conditions  on 
the  load-displacement  curve  and  the  buckling  load. 

e.  Effect  of  creep  on  the  buckling  load. 

Most  of  the  tests  were  performed  with  control  of  the  volume 
enclosed  under  the  shell  and  were  therefore  so-called  constant  volume 
tests .  This  name  comes  from  the  fact  that  the  enclosed  volume  does 
not  change  during  buckling.  Such  a  restriction  can  be  accomplished  by 
having  an  incompressible  fluid  filling  the  space  under  the  shell.  This 
type  of  test  was  chosen  so  that  a  load-displacemont  curve  could  be 
determined  since  both  pressure  and  a  measure  of  deflection  were  known 
at  any  time.  The  change  in  enclosed  volume  is  proportional  to  the 
change  in  average  deflection  of  the  shell  under  an  increment  of 
pressure.  The  constant  volume  test  was  also  chosen  so  that  the 
buckling  process  could  be  partially  controlled  and  forced  to  occur 
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gradually.  In  addition  it  seemed  desirable  to  have  either  constant 
pressure  or  constant  volume  conditions  so  that  the  results  could  be 
applied  to  the  discussion  of  the  energy  criterion  of  Tsien.  (in  a 
constant  pressure  test ,  the  pressure  remains  constant  during  buckling.) 

In  a  practical  situation  the  conditions  often  approximate  those 
of  a  constant  pressure  test.  Consider  a  hemispherical  structure  under 
a  wind  loading.  If  buckling  occurs  and  a  dimple  forms,  the  internal 
air  may  be  compressed  and  produce  some  small  resistance  to  volume  change. 
Thus  the  buckling  action,  although  mainly  a  constant  pressure  situation, 
is  really  somewhere  between  constant  pressure  and  constant  volume. 
Therefore  it  was  necessary  to  make  some  constant  pressure  tests.  The 
latter  tests  were  made  with  the  same  conditions  as  the  constant  volume 
tests  so  that  a  relation  between  buckling  loads  in  the  two  cases  might 
be  found.  With  such  a  relation,  buckling  loads  in  practical  situations 
could  be  predicted  from  the  results  of  constant  volu.Tie  tests. 

In  the  usual  structure  pressure  is  not  uniformly  distributed  over 
the  shell  but  may  have  a  variety  of  patterns.  However,  the  rate  of 
change  of  pressure  over  a  shell  surface  is  usually  small  and  buckling 
would  probably  occur  under  a  region  of  maximum  pressure.  If  only  the 
small  segment  which  forms  a  buckle  is  considered,  there  is  nearly  a 
uniform  pressure  applied.  Therefore,  it  is  reasonable  to  test  using 
a  uniform  pressure.  Also  the  uniform  case  is  easier  to  reproduce  and 
easier  to  compa.re  with  other  test  results.  A  uniform  pressure  was 
therefore  used  in  these  tests. 
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A  plastic  was  chosen  as  the  shell  material  rather  than  aluminum  or 


other  light  metals  which  have  usually  been  used.  There  were  a  number 
of  reasons  for  preferring  a  plastic  material.  A  plastic  was  relatively 
inexpensive  and  easy  to  form  by  vacuum  drawing;  it  had  reasonably  small 
thickness  variation;  and  there  was  no  strain-hardening  nor  residual 
stresses  developed  during  the  forming  operation.  It  was  desired  to  test 
entirely  in  the  elastic  range  of  the  material  so  that  each  shell  could 
be  buckled  many  times  and  not  cause  material  failure  or  yielding.  The 
large  ratio  of  yield  stress  to  modulus  of  elasticity  for  the  plastic 
allowed  the  large  deformations  of  buckling  to  occur  while  the  strains 
remained  elastic.  An  added  benefit  was  that  the  pressures  used  for  test 
ing  the  plastic  shells  were  l/20  to  l/60  of  that  for  comparable  metallic 
shells.  With  the  lower  pressures,  the  equipment  was  simpler,  less 
expensive,  and  more  convenient  to  use.  One  difficulty  with  the  plastic 
material  was  its  creep  characteristic.  The  nature  of  this  creep  was 
investigated  so  tha.,  test  results  could  be  reported  in  a  uniform  and 
meaningful  manner  and  extended  to  apply  to  ncn-creeping  materials.  The 
particular  plastic  chosen  was  a  polyvinyl  chloride  which  was  available 
in  thicknesses  from  0.010"  to  0.080". 

As  mentioned  in  Chapter  3,  the  form  of  deflection  seems  to  change 
at  certain  values  of  A  .  To  study  the  effect  of  the  mode  changes  on 
critical  loads,  it  was  decided  to  test  a  number  of  shells  with  A  values 
at  3.2,  5.3,  8.8,  16,  and  25,  that  is,  at  the  supposed  changeover  points 
As  a  check,  shells  with  A=  7.0,  12  and  20  were  also  used.  To  study 
reproducibility  of  loads  and  variations  of  loads  between  shells,  a  large 
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number  of  shells  were  formed  for  X=  7,0  and  12.  In  Fig.  13  it  is 
apparent  that  the  tests  are  grouped  about  certain  values  of  A .  That 
the  results  are  in  groups  instead  of  a  single  line  shows  that  the  desired 
^'s  were  only  achieved  approximately.  Procedures  described  below  were 
used  to  find  the  dimensions  of  the  shells  after  forming  and  these  dimen¬ 
sions  were  used  in  co.mputing  the  X '^or  each  shell.  After  choosing  X’s, 
it  was  necessary  to  select  values  of  thickness  and  radius  which  would 
give  the  desired  ^s.  This  choice  was  lim-* ‘-.ed  by  specifying  that  the 
buckling  pressure  must  be  greater  than  O.OSpsi.  and  that  the  sum  of 
bending  and  extensional  stress  at  the  edge  must  be  less  than  500  psi 
at  q  =  ^  q^.  From  Fig.  3  it  is  apparent  that  these  restrictions  left 
only  a  narrow  region  for  acceptable  combinations  of  h  and  a.  The 
possibility  was  considered  that  X  was  not  the  only  necessary  shape 
parameter.  Consequently,  where  possible  alternate  combinations  of  h 
and  a  were  made  for  a  single  A  .  The  small  circles  in  Fig.  3  show  the 
combinations  which  were  chosen.  For  each  set  several  shells  were  formed 
and  tested.  In  the  experiments  a  total  of  39  shells  were  used. 

Due  to  the  presence  of  creep  the  test  procedures  are  quite  important 
and  are  described  in  detail.  The  test  results  are  discussed  with  a  view 
to  answering  some  of  the  questions  which  have  been  raised  by  earlier 
tests  and  analyses  and  which  are  mentioned  in  Chapter  3. 
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B  Preparatory  Considerations 


-  Shell  Forming  - 

The  shell  segments  were  formed  from  polyvinyl  chloride  (P.V.C. ) 
sheets  by  a  process  known  as  vacuum  drawing.  In  this  method  a  single 
mold,  the  female,  was  required.  The  plastic  sheet  was  heated,  drawn  into 
the  mold,  and  allowed  to  cool  in  the  new  shape. 

Since  the  shells  were  to  have  five  different  radii  —  15",  25", 

35",  45",  and  80"  —  as  shown  in  Fig.  3,  five  molds  were  needed.  These 
molds  were  made  from  aluminum  plate  with  a  thickness  at  least  i"  greater 
than  the  height  of  the  shell  to  be  formed.  In  Fig.  4A  there  is  a  sketch 
cif  one.  such  mold.  Aluminum  is  relatively  easy  to  machine  but  has  better 
surface  and  is  stronger  than  aumy  of  the  other  possible  mold; materials 
such  as  wood  or  graphite.  A  cavity  with  one  of  the  desired  radii  was 
machined  i'nltoi  the  plate  and  a  l/l6"  diameter  hole  was  drilled  in  the  center 
of  the  cavity  for  access  to  the  vacuum. 

The  rest  of  the  forming  equipment  can  be  seen  in  Fig.  4A  also. 

The  wooden  clamping  rings  held  the  sheet  in  place  during  forming.  The 
large  cylinder  was  needed  to  provide  a  vacuum  chamber  and  to  hold  the 
sheet  slightly  above  the  mold  so  that  the  mold  was  Insulated  while  the 
sheet  was  being  heated.  The  cylinder  was  located  in  a  large  oven  and 
a  line  to  a  vacuum  pump  was  connected  to  control  pressure  in  the 
cylinder. 

The  forming  process  began  with  the  cutting  from  a  polyvinyl 
chloride  sheet  of  a  circular  piece  with  the  diameter  of  the  wooden 
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clamping  rings.  The  piece  was  rinsed  in  hot  running  water  to  remove 

grease,  dust,  and  other  foreign  particles.  Wiping  the  sheet  clean 

was  not  feasible  due  to  electrostatic  effects,  yet  cleaning  was 

necessary  or  bumps  would  appear  during  forming.  The  piece  was  clamped 

into  the  wooden  rings  with  negator  clamps  a:;d  placed  in  the  oven  atop 

the  cylinder  as  shown  in  the  lower  sketch  of  Fig.  4A.  With  the  sheet 

o  o 

to  be  formed  in  place,  the  oven  was  heated  to  150  -  170  C.  The  sheet 

seemed  to  be  properly  heated  when  the  steam  from  the  rinsing  water 
had  disappeared  from  it.  Then  the  vacuum  of  one  atmosphere  was  applied 
and  the  sheet  immediately  dropped  %o  fit  the  mold.  The  heat  was 
turned  off  and  the  oven  opened  to  speed  the  cooling.  About  one  minute 
was  allowed  for  anneailing  and  then  cold  water  was  poured  in  to  fill 
the  cavity  in  the  center  and  the  gutter  which  formed  just  outside  the 
mold  (see  the  second  sketch  in  Fig.  4A).  The  newly  formed  shell  was 
allowed  to  cool  under  full  vacuum  for  5  to  15  minutes.  The  longer 
period  of  cooling  was  required  for  the  sheets  0.080"  thick  while  the 
thinnest  sheets  needed  only  a  very  short  time.  After  the  cooling  process 
the  vacuum  was  turned  off  and  the  shell  was  lifted  out  of  the  mold. 

Then  shell  and  rings  were  dashed  together  into  a  vat  of  cool  water  so 
that  the  shell  could  cool  uniformly,  and  the  rings  were  removed.  At 
this  time  the  mold  was  also  taken  out  of  the  oven  so  that  it  would  be 
cool  for  the  next  forming. 

The  finished  shell  was  then  examined  to  find  any  visible  imperfections 
caused  by  the  presence  of  steam,  etc.  If  the  fault  was  in  the  forming, 
the  cycle  could  be  restarted  by  rinsing,  reclaraping  and  heating  as  before. 
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On  heating,  the  shell  became  planar  again  and  was  ready  for  redrawing. 

If  the  shell  passed  inspection,  the  excess  edge  material  was  trimmed 
off  so  that  only  3/4"  of  flat  rim  remained  and  the  specimen  was  ready 
for  testing  (see  Fig.  4A). 

The  vacuum  drawing  process  described  above  gradually  evolved 
during  the  production  of  the  shells  for  this  research.  Hence  some  of 
the  imperfections  in  the  first  shells  produced  were  eliminated  in  later 
production.  The  rim  slope  (see  Fig.  7)  was  caused  by  too  short  a  cooling 
time  under  vacuum  and/or  not  enough  heat.  Watermarks  (shallow  bumps) 
were  caused  by  trapped  steam.  Either  of  these  imperfections  could  have 
been  corrected  by  reheating  and  redrawing. 

-  Shell  Shape  Determination  - 

After  forming  the  shells  it  was  necessary  to  find  the  thickness 
and  radius  to  which  the  shell  actually  conformed.  Thicknesses  were 
read  to  the  nearest  ten  thousandth  of  an  inch  with  an  Ames  Dial  gage 
at  five  positions  in  the  shell  -  one  at  the  apex  of  the  shell,  and 
the  other  four  at  points  halfway  between  the  apex  and  the  edge.  An 
average  of  the  five  readings  was  taken  as  the  shell  thickness.  The 
thickness  variation  was  about  except  for  the  very  .thin  shells 
where  variations  were  10^  to  12^. 

The  radius  of  the  spherical  shell  can  be  found  if  the  rise  in  the 
center  is  known.  The  rise  was  measured  to  the  nearest  thousandth  of  an 
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inch,  Radius  was  then  computed  from  the  formula 


CL  ^ 


(4-1) 


where  H  is  the  central  rise  of  the  shell,  and  r  =  5.25”,  the  radius  of 
the  supported  circular  edge  as  seen  in  plan.  Unfortunately  some  of  the 
rims  of  the  shells  were  not  flat  but  had  a  slope  either  upward  or 
downward  as  shown  in  Fig.  7.  The  radius  determinations  for  these  shells 
were  nade  while  the  shells  were  clamped  into  the  testing  cylinders  in 
order  to  obtain  the  "a  tested"  radius  if  possible. 


-  Creep  Effect  on  Buckling  - 

The  creep  behavior  of  the  plastic  material  used  in  the  experiments 
is  similar  to  that  found  in  most  plastics.  The  total  creep  strain  at 
low  stress  levels  was  fairly  small,  around  5^  (at  room  temperature). 

The  creep  rate  was  very  high  so  that  creeping  was  essenti ally  complete 
in  twenty  minutes.  The  creep  became  nonlinear  for  stresses  over  a  few 
thousand  psi ,  that  is,  the  creep  rate  became  a  function  of  stress  level 
as  well  as  of  time.  The  creep  rate  varied  somewhat  with  temperature 
and  humidity.  At  moderate  stresses  the  creep  could  be  termed  elastic 
since  the  strain  was  completely  recovered  after  some  time.  In  fact, 
the  unloading  curve  matched  the  loading  curve. 

For  the  polyvinylchloride  material  used  in  the  experiments  the 
following  formula  was  found  to  describe  the  stress-strain  relation 
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adequately : 


£  =  InCti-t)]  (4-2) 

In  order  that  £ =  0  at  t  ~  0,  tg  was  defined  by 


/ 


This  makes  tg  extremely  small  so  that  it  can  be  neglected  for  most 
purposes  but  it  makes  the  evaluations  of  the  integrals  more  reasonable, 
g  is  a  creep  coefficient  and  is  the  modulus  at  one  minute  after 
loading. 

In  order  to  find  the  effect  of  creep  on  the  critical  load,  certain 
assumptions  were  made.  As  some  investigators  have  suggested,  buckling 
occurs  when  a  critical  strain  is  reached  at  some  point  in  the  structure 
and  does  not  depend  directly  on  the  stress  level.  In  the  following 
derivation  the  strain  was  assumed  to  be  a  linear  function  of  displacement, 
This  condition  is  reasonably  correct  to  the  point  of  buckling,  after 
which  nonlinearities  become  important. 

The  applied  load  may  not  be  constant  but  vary  with  time  7^,  In 
this  case  Eq.  (4-2)  is  still  applicable  but  the  stress  and  strain 
must  be  replaced  by  the  infinitesimal  quantities  d£  and  d <7"  .  A 
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(4-3) 


^  ^ a  /n  -7"  y  t.jj 

0(7'  ^ 

In  this  equation  7’  can  assume  any  value  from  0  to  t.  The  total 
strain  at  the  critical  point  in  the  shell  at  any  time  £  can  be  found 
by  integrating  the  preceding  equation  with  respect  to  7  from  7*0 
to  7  =  ^. 

set)  •  ^ 

If  the  load  is  applied  at  a  constant  rate,  then  O' "  where  T  is 

the  time  at  buckling. 


7^  In  7 1,)] 


/ 


When  t  =  T,  the  time  at  buckling,  ^(t)  =  critical 

stress  is 


where 


_ _ 

tt) 


A  normalized  stress  can  be  defined  as 


The  normalized  stress  is  a  function  only  of  Ej  and  neither  of 

which  is  a  function  of  the  duration  of  loading  T,  or  the  creep  coefficient 

g;  therefore,  (  (T  )„  is  not  a  function  of  T  or  g.  If  stress  is 
cr  N 
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proportional  to  laod,  then  normalized  buckling  pressures  can  be  similarly 
defined . 

"  K(r,3) 

P  is  a  nondimensional  form  of  buckling  pressure  defined  as  the  critical 
cr 

pressure,  divided  by  the  linear  buckling  load,  q^. 

In  equation  (4-4),  q^^  is  the  actual  pressure  measured  dur  ng  the 
test  and  (q  )  can  be  interpreted  as  the  pressure  for  a  hypothetical 
test  with' a  duration  of  2.728  minutes;  that  is,  for  K(T,g)  = 

1 .000'/'-;r— 1 .000  (since  to  can  be  neglected).  If  the  theory  is 
correct,  the  normalized  pressure,  like  the  normalized  stress,  is  not 
a  function  of  either  T  or  g.  Hence,  this  value  of  pressure  can  be 
compared  to  the  buckling  loads  for  noncreeping  materials  for  which  g 
IS  zero. 

Since  creep  adds  a  certain  amount  of  nonlinearity,  it  is  of 
interest  to  find  the  effect  of  creep  on  the  shape  of  the  load-deflection 
curve  up  to  the  point  of  buckling.  An  estimate  of  this  effect  was 
obtained  using  the  constants  mentioned  previously  and  a  T  of  twenty 
minutes.  Then  points  on  the  stress-strain  curve  were  obtained  for 
various  values  of  t  from  the  equation 

set)  ^  [i  !n 

A  linear  extrapolation  of  the  initial  portion  of  the  curve  was  made  to 
obtain  a  stress  at  the  abscissa  of  critical  strain.  This  value  from  the 
linear  extrapolation  was  then  compared  with  the  critical  stress  obtained  from 


^(T.3) 
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the  equation.  The  difference  in  this  case  was  only  0,3^.  Therefore 
it  is  fairly  safe  to  say  that  creep  did  not  account  for  nonlinearities 
found  in  the  load-deflection  curves. 

-  Modulus  Determination  - 

In  order  to  relate  these  test  results  to  other  materials,  some 
procedure  had  to  be  found  for  evaluating  the  modulus  of  elasticity. 

Simple  tensile  tests  were  run  on  samples  of  each  sheet  using  SR-4 
strain  gages  and  a  dead  weight.  Beam  tests  were  also  run  in  an  attempt 
to  determine  the  compressive  modulus  in  case  it  differed  from  the  tensile. 
Both  sets  of  tests  were  conducted  on  6"  by  ij"  specimens.  In  the  beam 
tests  the  specimens  were  supported  as  simple  beams  with  a  load  at  the 
midpoints  and  the  end  rotations  were  measured  by  the  motion  of  a  wire 
fixed  to  the  beam  near  the  end.  Although  the  beam  deflections  were 
extremely  large ,  the  relation  between  load  and  end  rotation  was  essentially 
linear.  In  Appendix  B  it  is  shown  that  the  end  rotation  equals 


where  ^  is  the  end  rotation  in  the  linear  theory.  Hence  for  a  rotation 
of  6°  or  0.1  radians,  which  was  about  the  maximum  used,  the  nonlinear 
correction  was  less  than  0,2^.  Because  of  the  presence  of  creep,  the 
modulus  changed  with  time.  It  was  decided  to  use  Ej ,  the  modulus  at 
one  minute  after  loading.  is  defined  by  the  creep  equation 
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(T 

Since  to  is  negligible,  t  =  1  minute  gives  = -^  and  there  is  no  creep 
correction  required  for  the  determination  of  the  modulus.  The  beam 
tests  gave  the  most  uniform  results  for  modulus,  the  average  value  being 
498,000  psi .  A  summary  of  these  determinations  is  in  Appendix  B. 

The  tensile  tests  were  not  so  reliable  due  to  the  presence  of 
the  strain  gages.  On  such  small  specimens  the  stiffness  of  the  gage 
produced  an  important  effect  on  the  strains  and  hence  on  the  modulus 
of  elasticity.  However  the  electric  strain  gages  could  be  read  with 
more  precision  than  the  beam  rotations.  Therefore  the  gages  were 
valuable  in  determining  creep  which  is  a  small  variation  on  the  total 
strain.  Hence  the  creep  equation  (4-2)  and  creep  coefficient,  g,  are 
based  on  the  tensile  tests  mainly.  Results  from  these  creep  and 
modulus  tests  are  also  given  in  Appendix  B. 
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C  Testing 


-  Test  Apparatus  - 

The  main  tests  were  constant  volume  tests  which  required  the  presence 
of  an  incompiessible  fluid  in  a  chamber  below  the  shell  (shell  convex 
up).  Water  was  selected  as  the  most  ccnvenient  fluid  for  this  purpose. 

But  with  water  below  and  air  above  the  net  pressure  on  the  shell  would 
be  hydrostatic  and  vary  over  the  shell  surface.  Therefore  there  had  to 
be  a  water  chamber  above  the  shell  also.  Tests  were  run  by  applying 
a  displacement  and  measuring  the  pressure  taken  by  the  shell.  Hence 
the  required  apparatus  was  a  set  of  water  tight  chambers  between  which 
the  shell  could  be  clamped,  a  device  for  withdrawing  water  from  the 
lower  chamb»?r,  and  a  syste.m  for  measuring  pressures  above  and  below  the 
shell . 

The  completed  equipment  can  be  seen  in  Fig.  4.  The  lower  chamber 
consisted  of  a  heavy  pipe  section  bolted  to  a  base  plate.  These  two 
pieces  were  rather  heavy  so  that  they  would  not  alter  the  enclosed 
volume  much  as  the  pressures  changed.  The  upper  chamber  was  made  of 
plexiglas  and  allowed  observation  of  the  shell,  the  buckles,  and  bub¬ 
bles  which  were  trapped  in  the  chambers.  Holes  for  inlets  and  out¬ 
lets  were  made  at  three  places  around  the  lower  chamber.  These  were 
necessary  in  crder  to  allow  the  water  to  come  in  and  the  air  bubbles 
to  be  forced  out,  It  was  found  convenient  to  use  water  from  a  large 
jug  which  maintained  a  nearly  constant  water  temperature  so  that  the 
volume  would  not  change  due  to  temperature  fluctuations. 
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The  tubing  which  connected  the  reservoirs,  chambers,  and  mano¬ 
meters  was  of  transparent  plastic  which  allowed  observation  of  bubbles 
in  the  line.  When  bubbles  were  present, the  manometer  readings  were 
very  unsteady  since  each  motion  of  the  bubbles  seemed  to  upset  the 
equilibrium* cf  the  gages,  too. 

The  manometers  were  capillary  tubes  with  an  inside  diameter  of 
1.8  mm.  These  were  used  so  that  when  the  pressure  changed,  the  change 
in  the  amount  of  liquid  in  t)ie  tube  would  not  greatly  alter  the  total 
volume  under  the  shell.  The  tubes  for  registering  upper  and  lower 
pressure  had  equal  diameters  so  that  the  capillary  rise  (i"  to  l") 
would  be  the  same  and  no  correction  was  required.  The  capillary 
rise  was  quite  a  bit  larger  in  pure  water  than  in  impure  and  so  small 
amounts  of  liquid  soap  were  added  in  both  tubes  to  make  certain  the 
water  was  impure.  The  soap  reduced  wall  friction  between  glass  and 
water  and  so  also  produced  freer  and  more  accurate  motion  of  the 
water . 

The  faucet  for  allowing  withdrawal  of  water  from  the  lower  chamber 
was  at  first  a  high  vacuum  glass  valve.  Later  tbs  valve  was  replaced 
by  a  short  rubber  tube  with  a  small  clamp  on  it.  With  the  clamp  there 
was  good  control  of  the  rate  at  which  water  drops  came  out  so  that  a 
desired  speed  of  volume  change  could  be  maintained. 

The  rubber  gaskets  used  in  clamping  the  shell  may  have  caused 
some  inaccuracies.  The  narrow  rim  all  around  each  shell  was  supposed 
to  be  clamped  against  all  motion.  However,  the  thick  gaskets  in  the 
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clamping  rings  probably  allowed  some  outward  displacement  and  rotation 
of  the  edges.  Rough  calculations  show  that  the  amount  of  displacement 
would  be  quite  small  compared  to  the  total  shortening  of  a  meridian  and 
is  therefore  not  too  important  as  a  cause  for  snapping  through.  The 
edge  rotation  would  not  approach  that  for  the  simply  supported  case, 
yet  the  theoretical  buckling  load  for  the  latter  c?.se  is  at  least  3/4ths 

of  that  for  the  fixed  condition  according  to  the  calculations  of 

^  (12b) 

Weinitschke 


-  Test  Procedure  - 

There  were  two  basically  different  types  of  tests  used.  The 
con.stant  volume  test  was  the  most  important  and  it  was  of  the  con¬ 
trolled  displacement  type.  It  proceeded  much  like  an  ordinary  tensile 
test.  A  certain  strain  or  displacement  was  applied  to  the  shell  and 
the  pressure  which  was  required  to  make  the  shell  stay  in  that  position 
was  read.  The  other  type,  the  constant  pressure  test,  was  carried  out 
by  increasing  the  load  until  buckling  occurred  and  hence  the  control 
was  on  the  load,  not  on  the  displacement. 

Preparation  for  the  constant  volume  tests  was  commenced  by  placing 
the  shell  betwee.’?  the  clamping  rings  of  the  test  cylinders  and  bolting 
them  together,  as  shown  in  Fig.  4.  Then  the  chambers  were  filled  with 
water,  beginning  with  the  lower  chamber.  The  tube  to  the  top  chamber 
was  connected  to  the  small  elevated  reservoir  and  the  tube  to  the  lower 
chamber  was  connected  to  the  taller  capillary  tube  as  shown  in  Fig.  4. 
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Equal  heights  were  obtained  in  the  upper  and  lower  chamber  manometers 
by  letting  water  either  in  or  out  of  the  lower  chamber  through  the  faucet. 
To  force  water  in,  the  reversing  standpipe  shown  in  Fig.  4  was  attached 
to  the  valve  and  fil'ed  with  water.  With  equal  pressure  above  and  below 
the  shell,  the  shell  was  sustaining  a  net  load  of  zero. 

Testing  began  by  opening  the  faucet  and  allowing  a  number  of 
drops  of  water  to  come  out  from  the  lower  chamber.  The  shell  had  to 
move  downward  in  order  to  keep  in  contact  with  the  water  and  thus  the 
displacement  was  applied  to  the  shell.  Of  course,  a  certain  pressure 
was  required  to  hold  the  shell  in  the  deflected  shape  and  this  pres¬ 
sure  was  given  by  the  difference  between  pressures  in  the  upper  and 
lower  chambers.  Hence  the  load  on  the  shell  could  be  found  from  the 
readings  of  the  two  tManometers .  So  the  normal  procedure  was  to  open 
the  tap,  let  out  a  number  of  drops,  close  the  tap,  read  manometers, 
and  then  repeat.  The  information  recorded  in  each  cycle  was  the 
number  of  drops  and  both  manometer  readings.  At  buckling  the  height 
in  the  lower  chamber  manometer  increased  suddenly  and  readi.ngs  were 
taken  both  before  and  after  the  Jump.  Then  the  usual  procedure  was 
continued  on  to  some  point  in  the  post-buckling  range.  For  deep 
shells  the  qa ,  the  pressure  to  which  the  shell  jumped  at  buckling, 

and  q  ,  ,  the  pressure  at  the  minimum  point  of  the  load-deflection 
min 

curve  (see  Fig,  2),  coincided  and  so  the  test  was  often  halted  just 

after  buckling.  For  shallower  shells,  the  curve  looked  more  like  the 

one  in  Figure  2  and  the  test  was  continued  until  q  had  been  passed. 

min 

See  Figure  10  for  some  typical  load-displacement  curves. 
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After  the  test  the  total  amount  of  waiter  taken  out  during  the 


test  ^as  determined  by  weighing.  This  information  was  used  to  associate 
a  volume  change  with  each  pressure  reading  by  apportioning  the  total 
change  according  to  the  total  drops  up  to  each  point.  For  plotting 
load  versus  volume  change  curves,  it  was  necessary  only  to  plot  the 
pressure  versus  the  cumulative  number  of  drops.  The  drops  were  found 
to  be  somewhat  nonuniform  in  size  but  not  enough  to  make  the  initial 
portion  of  the  load-deflection  curve  appear  nonlinear.  Therefore, 
the  number  of  drops  was  probably  an  adequately  reliable  measure  of 
volume  change. 

In  order  to  unbuckle  the  shell  the  reversing  column  was  attached 
to  the  faucet.  Then  the  normal  testing  procedure  was  to  open  the  tap, 
allow  a  certain  amount  of  water  to  run  into  the  lower  chamber,  close 
the  tap,  and  read  the  heights  of  manometers  and  reversing  column.  The 
difference  between  the  current  reversing  column  reading  and  the  previous 
one  gave  the  volume  change  because  the  column  was  calibrated.  The 
procedure  was  repeated  until  some  point  on  the  initial  portion  of  the 
load-def lectio.n  curve  was  reached. 

With  the  presence  of  creep,  time  became  a  factor  in  the  test. 

The  steps  of  the  con.stant  volume  test  enumerated  below  were  adopted  to 
standardize  the  effect  of  time. 

a.  The  pressures  in  the  water  manometers  and  the  time  were 
recorded . 

b.  The  tap  was  opened  to  allow  a  number  of  drops  (say  10)  to 
escape  fiom  the  lower  chamber  and  rhe  number  was  recorded. 
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Again  the  pressures  in  the  gages  were  read  and  noted.  At 


c  . 

.some  specified  time  (say  30  seconds)  after  the  last  time  the  tap  was 
opened,  step  b  was  repeated.  By  using  an  equal  time  interval  during 
each  cycle,  the  loading  rate  was  made  to  be  constant  stepwise. 

d.  At  buckling  the  tap  was  immediately  turned  off  and  the  time, 
pressure  readings,  and  number  of  drops  in  the  last  cycle  were  recorded. 
The  remainder  of  the  experiment  proceeded  without  being  timed. 

e.  Steps  b  and  c  were  repeated  until  the  desired  point  on  the 
post-buckling  curve  had  been  reached. 

The  loading  process  was  then  reversed  by  attaching  the  reversing 
column  to  the  tap  and  filling  it  with  water.  In  the  rebound  test  a 
procedure  similar  to  that  above  was  followed. 

a.  The  pressures  in  the  water  manometers  and  the  height  of  the 
reversing  column  were  read  and  recorded. 

b.  The  tap  was  opened  to  allow  a  number  of  centimeters  of  water 
trom  the  column  to  enter  the  lower  chamber. 

c.  Steps  a  and  b  were  repeated  until  the  shell  had  unbuckled 
aiid  the  pressure  was  reducing  along  the  initial  portion  of  the  load- 
deflection  curve.  See  Fig.  10  for  some  typical  rebound  curves. 

The  preceding  tests  give  a  complete  curve  of  load  versus  displace- 
.meni  measured  in  terms  of  the  change  in  volume  under  the  shell.  After 
the  nature  of  this  curve  was  known  for  a  shell  or  a  group  of  similar 
shells,  the  only  things  that  were  required  of  further  tests  were  certain 
significant  pressures  on  the  curve  such  as  q  ,  q2>  Q  .  and  q  .  .  The 
following  modified  procedure  gave  these  impc- tant  pressures  and  was  used 
for  a  majority  of  the  tests. 
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The  time  was  noted  and  the  tap  was  opened, 

b.  A  uniform  drop  rate  was  maintained.  This  was  done  visually 
for  short  buckling  time^  but  for  long  tests,  the  rate  was  kept  oven 
by  passing  an  equal  number  of  drops  each  minute. 

c.  While  keeping  a  uniform  drop  rate  the  pressure  in  the  lower 
chamber  manometer  was  observed  continuously. 

d.  When  buckling  occurred,  the  faucet  was  shut  off,  and  the 
time  and  the  reading  of  the  lower  manometer  prior  to  buckling  were 
recorded.  The  pressure  in  the  upper  chaunber  manometer  was  read  and 
recorded.  At  one  minute  after  buckling  the  height  of  the  lower 
manometer  was  read  again.  This  latter  reading  was  taken  to  find  the 
so  called  lower  buckling  pressure,  q2 t  Just  after  buckling  the  height 
in  the  lower  manometer  was  changing  rather  rapidly  due  to  creep  so  the 
specification  of  a  one  minute  wait  was  necessary  to  make  the  readings 
uniform. 

e.  The  faucet  was  turned  on  again  and  the  height  of  the  lower 

manometer  was  observed.  If  the  height  increased  and  then  decreased 

(pressure  on  the  shell  decreased  and  then  increased)  then  a  reading 

was  taken  at  the  maximum  to  give  q  ^  .  Otherwise  it  was  assumed  that 

min 

qa  =  q  .  . 

^  min 

f.  The  reversing  column  was  used  as  before  to  force  water  back 
into  the  lower  chamber.  The  tap  was  turned  on  and  the  lower  chamber 
manometer  was  observed. 
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The  warning  for  the  onset  of  unbuckling  was  a  gradually 


accelerated  decrease  in  gage  height.  Near  the  time  of  unbuckling  the 
tap  was  turned  off  and  the  shell  was  allowed  to  unbuckle  without 
outside  influence. 

h.  At  one  minute  after  unbuckling,  the  two  gage  readings  were 
taken  and  recorded.  The  creep  in  gage  height  was  generally  not  large 
in  this  case  but  still  noticeable. 

After  the  pressures  in  the  two  chambers  were  equalized  following 
the  test,  a  period  of  20  minutes  was  allowed  to  elapse  before  the  next 
test  was  commenced.  This  wait  time  seemed  necessary  so  that  the  shell 
would  have  essentially  no  memory  (due  to  creep)  of  the  previous  test. 

In  order  to  check  the  reproducibility  of  buckling  loads,  a  series 
of  tests  was  run  without  removing  the  shell  from  the  testing  cylinders. 
Then  the  shell  was  taken  out  and  replaced  in  the  test  chambers  and 
tested  again.  This  procedure  also  served  to  demonstrate  the  effect 
of  precision  in  clamping  the  shell  in  position. 

After  a  series  of  constant  volume  tests  had  been  run  on  some  of 
the  shells,  the  water  was  emptied  out  and  a  "constant  pressure"  test 
was  run.  Actually  the  volume  of  air  was  not  great  enough  to  allow 
for  full  collapse  of  the  shell  under  the  same  pressure,  but  the  test 
approximated  constant  pressure  conditions  to  some  degree.  This  degree 
may  be  estimated  from  the  fact  that  a  lower  buckling  load  was  obtained 
and  that  the  buckle  was  about  the  same  size  as  that  for  constant  volume 
tests.  The  maiii  distinction  then  was  that  the  pressure  was  applied. 
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not  the  displacement.  In  these  tests  a  U  tube  manometer  and  pressure 
inlet  were  both  attached  to  the  upper  chamber.  The  lov^er  chamber  was 
open  to  the  atmosphere.  Pressure  was  built  up  gradually  and  the 
manometer  was  watched  continuously.  When  buckling  occurred,  the 
maximum  manometer  height  and  duration  of  the  test  were  recorded. 
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D  Experimental  Results 


-  General  Discussion  - 

The  data  from  both  constant  volume  and  constant  pressure  tests 

are  given  in  Table  1  of  Appendix  A.  The  tests  reported  were  made 

in  accordance  with  the  procedures  just  described.  Individual  test-s 

are  not  reported,  only  series  of  tests  which  were  made  under  the  same 

conditions.  For  instance  "a"  is  the  designation  for  a  series  of 

constant  volume  tests  and  "a(air)"  is  a  series  of  constant  prefisure 

tests  run  under  the  same  clamping  conditions  as  the  constant  volume 

tests.  The  values  of  P  in  each  series  were  normalized  according 

cr 

to  Eq.  (4-4)  and  then  averaged.  This  procedure  removed  the  effect 
of  different  lengths  of  time  in  testing.  Under  series  "a"  were 
reported  the  averaged  values  of  (P  )  ,  Po ,  P  and  the  position  of 
the  buckle  whicn  formed.  Notice  that  during  any  series  the  buckle 
occurred  in  the  same  position.  The  one  exception  to  this  is  the 
second  series  of  tests  on  number  35. 

The  buckle  and  imperfection  positions  were  given  as: 
c  center  of  shell 

n  near  center 

m  midway  between  center  and  edge 

e  adjacent  to  edge 

For  each  shell  the  values  of  q,^,,  ^  ,  a/h,  and  0^  are  listed.  This  is 

followed  by  a  brief  description  of  observed  imperfections.  First  was 
ihe  slope  of  the  rim  of  the  shell.  If  the  edge  strip  curled  up  (with 
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the  shell  convex  up)  after  forming  a  positive  s]ope  was  indicated. 

Next  was  the  estimate  of  the  height  of  bumps  and  their  location  on 
the  shell. 

From  any  series  of  tests  run  under  the  same  conditions,  it  was 

found  that  the  value  of  q  was  nearly  the  same  in  all  tests.  When 

cr 

the  buckling  pressures  were  normalized  according  to  Eq .  (4-4)  to 
account  for  creep,  then  the  variations  were  less  than  w' thin  any 
one  series.  With  such  small  variations  it  was  possible  to  conclude 
that  buckling  pressures  were  reproducible  and  dependent  on  conditions 
which  were  actually  unvaried  during  a  series  of  tests.  Between  series 
of  tests  there  was  much  more  spread  in  the  buckling  loads,  but  the 
variations  were  still  of  the  order  of  5  or  10^,  which  is  quite 
reasonable.  The  greatest  variations  occurred  between  different  shells 
with  supposedly  the  same  shape.  In  a  number  of  cases  the  critical 
load  for  one  shell  was  50^  of  that  for  another  with  the  identical  gross 
dimensions.  A  number  of  test  series  were  made  with  large  variations 
in  test  durations  in  order  to  check  on  the  creep  correction,  Eq.  (4-4). 
One  such  series  is  plotted  in  Fig.  6  and  from  the  correlation  shown 
there  between  theory  and  experiment,  it  is  apparent  that  the  correction 
is  of  about  the  right  magnitude. 

The  buckling  pressures  from  the  air  pressure  tests  averaged, 
about  3^  higher  than  those  from  constant  volume  tests.  This  comparison 
is  made  between  results  from  experiments  all  made  without  changing 
clamping  conditions.  That  the  conditions  were  quite  constant  between 
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the  two  types  of  tests  is  shown  by  the  fact  that  the  buckle  formed  in 
the  same  position  in  both  cases.  No  explanation  has  been  found  for  the 
increase  in  load  under  constant  pressure  conditions. 

The  value  of  A  for  the  shells  were  chosen  at  the  points  where 
the  modes  of  deflected  shape  were  expected  to  change,  and  at  intermediate 
points.  The  buckle  position  does  not  seem  to  depend  noticeably  on 
whether  or  not  ^  is  at  a  changeover  point  or  not  (Fig.  9a).  Neither 
IS  there  any  particular  difference  in  buckling  pressures  between  shells 
in  the  two  groups. 

The  nature  of  qj ,  the  pressure  to  which  the  shell  jumped  at  buckling, 
and  the  pressure  found  at  unbuckling,  contrasted  markedly  with  that 

of  q  .  (qo  and  q  are  from  constant  volume  tests  only.)  The  spread 
between  qa  and  was  not  very  great  as  shown  in  Fig.  13.  In  fact  it 
was  reasonable  to  draw  in  the  trend  lines  shown.  In  view  of  this  uniformity 
of  averaged  results,  it  was  surprising  to  find  that  qg  and  q^^  were 
actually  not  very  reproducible  under  the  same  conditions.  For  instance, 
in  a  single  test  series  in  which  there  was  little  variation  in  q  , 
the  spread  in  values  of  q^  might  be  10%  or  even  30%.  The  variations 
between  values  of  qa  in  different  series  was  not  greater  than  the 
variation  within  the  series.  Apparently  all  the  conditions  which 
affect  the  value  of  qo  and  q  were  not  under  control.  No  corrections 
were  made  for  creep  and  that  may  have  been  a  partial  cause  of  the 
variations . 
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For  the  deep  shells  there  was  a  long,  level  portion  in  the  post- 

buckling  curve  and  q,  and  q  both  fell  on  this  part,  so  that 

min 

q-  =  q  .In  the  shallower  shells  qa  was  sometimes  quite  different 
min  * 

from  q  ,  .  The  upper  sketch  in  Fig.  10  is  for  such  a  shallow  shell, 
min 

Not  all  of  the  important  pressures  --  *^min 

depended  solely  on  A  ,  although  this  shape  parameter  was  dominant  in 

most  cases.  From  the  plot  of  P  versus  X  in  Fig.  13a  it  can  be  seen 

that  P  =  q  /qT  was  essentially  independent  of  A  and  there  was  enough 
cr  Clr 

spread  in  the  data  to  obscure  any  minor  trends.  There  is  just  a  slight 

tendency  for  P  to  decrease  with  increasing  X  .  q,  and  q  seemed  to 
cr  *  un 

be  related  more  closely  to  A  .  as  evidenced  by  the  trends  seen  in 

Fig.  13.  In  fact,  q2  and  q^^  were  replotted  versus  0('  and  versus  a/h, 

but  much  better  correlation  was  found  in  the  plots  with  X  as  abscissa. 

Since  many  analysts  have  predicted  a  negative  value  for  (which 

equals  qj  for  large  A  )i  it  might  be  expected  that  qj  would  at  lea.st 

tend  to  zero  as  X  increases.  Extrapolation  is  risky  with  such  scattered 

results  but  the  trend  seems  to  i.ndicate  that  the  curve  of  q^  would 

become  asymptotic  to  zero  but  not  reach  a  negative  value. 

The  value  of  q  (from  constant  volume  tests)  seemed  to  be 
min 

related  to  sh"ll  height  instead  of  to  X .  This  relationship  was 
indicated  by  the  differences  in  behavior  between  shells  8  and  14. 

These  two  have  about  the  same  value  of  A  ,  but  the  load-displacement 
curves  associated  with  them  are  quite  different.  In  Fig.  10  it  is  seen 
that  after  buckling,  the  pressure  continued  to  decrease  as  the  deflection 
increased.  Thus  '*'as  much  less  than  q2 .  But  v  j  th  shell  14,  the 
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curve  st&rted  up  again  immediately  after  buckling  so  that  Qa  = 

Probably  q  should  be  plotted  with  o<  as  thu  abscissa  for  good 
min 

correlation.  Notwithstanding  what  has  just  been  said,  it  appeared  from 

Fjg.  9b  that  P  -  P  .  is  probably  a  function  of  A.  The  same  plot 
cr  min 

was  attempted  with  P^^  -  Pj  and  the  scatter  was  worse  in  the  low 

p:!itlon  which  determines  the  relationship.  However,  the  conclusion 

that  ?  -  P  IS  a  functionof  X  >nay  have  been  due  to  a  fortuitous 

cr  min  j 

scatter  in  values  of  P 

cr 


-  Buckling  Process  - 

The  entire  load-displacement  curve  varies  a  great  deal  depending 
on  the  value  of  the  parameter  X  and  on  the  shell  depth.  For  the  very 
shallow  shells  (small  X)  the  curve  became  extremely  nonlinear  before 
buckling  and  was  essentially  horizontal  at  the  time  of  buckling.  A 
curve  for  a  fairly  shallow  shell  is  shown  in  the  upper  sketch  of  Pig.  10. 
Buckling  then  occurred  by  a  gradual  decrease  in  the  pressure  maintained 
by  the  shell.  The  pressure  after  buckling  was  only  a  few  percent  below 
tht^  maximu.m  or  buckling  pressure  of  the  shell.  For  the  case  of  large 
X  tr.e  situation  was  entirely  different.  Here  the  load-displacement 
car.e  was  t.-&sentially  linear  up  to  the  point  of  buckling.  Buckling  was 
extremely  sudden,  and  without  warning.  The  post-buckling  pressure  was 
about  if-i  to  1/6  of  the  maximum.  A  curve  for  a  deep  shell  is  also 
shown  in  Fig  10. 
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In  all  cases  buckling  was  accompanied  by  the  formation  of  one 
buckle  or  dimple  (see  Fig.  5).  The  buckle  under  water  pressure  was 
smooth,  not  jagged  as  in  the  constant  pressure  tests  of  Kloppel  and 
Jungbluth'  The  total  central  angle  of  the  buckle  is  about  7  for 

all  X.  Actually  the  size  of  the  buckle  is  a  peculiarity  of  the  constant 
volume  test  and  not  of  the  buckling  process  or  of  the  shell.  When 
further  volume  was  withdrawn  after  the  occurrence  of  buckling,  the  size 
of  the  buckle  increased.  When  the  buckle  was  at  maximum  size,  it 
extended  over  the  entire  shell.  When  still  more  volume  was  taken  out 
after  the  buckle  reached  maximum  size,  the  pressure  taken  by  the  shell 
began  to  increase.  If  volume  was  then  forced  back  into  the  shell,  the 
buckle  decreased  until  its  central  angle  was  about  2°  or  3*^  before 
snapping  back  through. 

The  buckle  could  occur  at  any  location  in  the  shell.  However, 
if  the  shell  was  not  removed  from  the  cylinders  between  tests  (no 
change  in  test  conditions),  then  the  buckle  occurred  again  in  the 
same  position  (e.xcept  in  test  series  35b).  But,  under  changed  boundary 
conditions,  as  when  the  shell  was  undamped  and  replaced  in  the  test 
apparatus,  it  was  not  possible  to  predict  the  new  location  of  the 
buckle . 

For  identification  purposes,  four  buckle  positions  were  chosen: 
near  the  edge,  midway  between  the  center  and  edge,  near  the  center, 
and  center.  By  correlating  buckle  positions  with  buckling  loads  for 
the  shells,  it  was  found  that  edge  buckles  were  usually  associated 
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with  low  loads,  center  buckles  with  high  loads,  and  the  other  two  have 
intermediate  buckling  loads.  This  correlation  seems  reasonable  from 
the  following  viewpoint.  In  any  given  shell  there  will  be  Some  imper¬ 
fections  and  Some  moments  induced  at  the  boundary  by  the  clamping  rings. 
It  the  effect  of  the  imperfection  dominates,  the  same  load  should  be 
obtained  in  all  tests,  end  the  buckle  should  be  formed  in  the  vicinity 
of  the  imperfection.  However,  if  the  edge  moments  are  domlttant,  the 
buckle  will  be  formed  near  the  edge  and  a  lower  load  must  be  found, 
Since  the  imperfections  will  not  be  changed  from  test  to  test,  the 
effect  of  any  other  influence  must  be  either  none  or  such  as  to  de¬ 
crease  the  buckling  load.  However,  in  this  reasoning  a  very  simple 
relation  between  imperfection  and  buckling  load  was  assumed. 

Unbuckling  was  accomplished  by  forcing  water  back  under  the  shell. 
For  the  deeper  shells  an  unbuckling  action  occurred  much  like  the  buck¬ 
ling.  At  some  point  there  was  a  sudden  increase  in  the  pressure  main¬ 
tained  by  the  shell.  This  occurred  as  the  buckle  snapped  back  through. 
Then  the  shell  was  back  on  the  initial  portion  of  the  load-displacement 
curve  and  further  volume  change  caused  a  decrease  in  the  pressure  in  a 
linear  manner.  Just  prior  to  unbuckling,  it  was  observed  that  the 
dimple  was  vibrating  at  three  or  four  cycles  per  second.  No  such 
motion  was  noticed  at  any  other  time  during  the  loading  or  unloading 
process . 

Various  investigators  have  determined  the  minimum  value  of  \ 
at  which  buckling  occurs.  Below  this  value  the  load  deflection  curve 
is  a  .single-valued  function  of  pressure,  Kaplan  and  Fung  found 
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2,1  theoretically,  but  5.0  experimentally  as  the  minimum  for 
buckling.  In  his  analysis  Weinitschke ^  determined  3.2  as  the 
limiting  value.  In  the  tests  of  this  report,  one  shell  with  X  =  3.3 
exhibited  an  inflection  point  but  no  decrease  in  load  with  increase  in 
deflection.  Another  shell  with  X=  3-4  showed  a  slight  decrease  in 
pressure  in  some  tests  and  none  in  others.  These  tests  are  therefore 
some  verification  for  the  theoretical  value  of  Weinitschke. 

-  Effects  of  Imperfections  - 

One  type  of  imperfection  encountered  in  the  shells  was  a  slope 
of  the  rim  as  depicted  in  Fig.  7.  The  sign  convention  adopted  was 
also  shown  in  this  figure.  Some  change  in  buckling  loads  might  be 
expected  due  to  this  slope.  When  the  shells  were  clamped  into  the 
test  chambers,  the  rims  were  forced  fist,  thus  inducing  moments 
along  the  shell  edge  and,  for  large  slopes,  the  shell  height  was 
even  altered.  Fig.  7  was  plotted  using  the  simple  assumption  that 
a  positive  slope  would  tend  to  increase  the  buckling  load  and  a 
negative  slope  would  have  tlie  opposite  ef-^ect.  However  no  trends 
at  all  are  indicated  in  the  figure. 

Hudiansky  made  an  observation  concerning  the  interaction  of 
mode  shapes  and  boundary  conditions  which  may  be  used  here.  During 
loading  a  trough  forms  about  the  shell  edge  as  a  transition  from  the 
clamped  boundary  to  che  uniform  deflection  which  occurs  over  the 
center  of  the  shell.  The  mode  shapes  are  also  a  series  of  troughs 
and  crests,  and  the  number  and  position  of  these  waves  depend  on  the 
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value  of  A  .  for  some  values  of  A  a  trough  should  appear  at  the  edge, 

for  others,  a  crest.  Budiansky  showed  that  when  the  mode  shape  produced 

a  trough  at  the  boundary,  reinforcement  occurred  with  the  clamping 

trough  and  low  loads  were  obtained  at  between  A  =  3.5  and  5.5  and 

between  8  and  12.  For  A  between  5.5  and  8,  cancellation  occurred 

and  buckling  pressure  was  higher.  As  A  increases  the  oscillation 

in  the  P  versus  A  curve  dies  out  because  the  conditions  at  the 
cr 

apex  and  edge  of  the  shell  have  less  mutual  influence.  Fig.  7  was 
replotted  using  these  ideas.  For  shells  with  A  between  3.5  and  5.5, 

8  and  12,  and  between  16  and  25  the  same  sign  convention  was  used  for 
the  rim  slope.  For  the  other  shells  the  signs  of  the  slopes  were 
reversed.  The  figure  which  was  thus  produced  looked  no  different  from 
Fig.  7.  That  no  trends  appeared  in  either  figure  should  not  be  taken 
as  a  sign  that  there  was  no  efect  from  the  slope  of  the  rim.  Rather 
it  was  an  indication  that  the  influence  of  slope  on  buckling  pressure 
is  small  and  obscured  by  other  things. 

Interest  in  the  effect  of  clamping  a  sloping  rim  is  not  actually 
restricted  to  this  set  of  tests.  In  a  practical  structure,  edge 
restraint  is  provided  generally  by  other  deformaoie  bodies  which  will 
tend  to  curl  the  shell  edge  up  or  down  in  a  manner  similar  to  the  above 
effect 

Various  surface  imperfections  were  found  on  the  shells.  Among 
these  blemishes  were  central  bumps  which  were  caused  by  an  irregu¬ 
larity  in  one  of  the  molds  and  miscellaneous  shallow  bumps  caused 


56 


by  the  presence  of  steam  under  the  shell  during  forming.  In  a  few 
shells  there  were  very  small  sharp  bumps  caused  by  dust  flacks  present 
during  forming,  A  further  imperfection  was  caused  if  the  hot  sheet 
did  not  hit  the  mold  uniformly  when  the  vacuum  was  applied  so  that  air 
pockets  remained.  Such  an  imperfection  was  a  symmetric  groove  at 
midheight  around  shell  16,  Shell  25  (and  to  a  lesser  extent,  24) 
had  a  large  fTat  section  at  the  apex.  This  seems  to  have  reduced  the 
load  considerably  and  even  flattened  the  load-deflection  curve  so  that 
buckling  occurred  gradually,  not  suddenly. 

Fig.  8  was  a  plot  of  Imperfection  magnitude  versus  the  buckling 
pressure,  No  trend  was  evident  but  again  this  was  based  on  a  very 
simple  idea  about  the  effect  of  imperfections.  It  was  not  felt  that 
a  more  valuable  plot  could  be  made  using  the  knowledge  which  is  presently 
at  hand.  The  analysis  of  Chen  has  shown  that  buckling  pressures  may  be 
either  increased  or  decreased  by  imperfections  defending  on  the  value 
of  ^  and  on  the  position,  size,  and  direction  of  the  imperfection.  This 
is  reasonable  if  one  considers  the  fact  that  shell  structures  are  usually 
built  with  ridges,  grooves,  and  flat  spots  for  added  rigidity  and 
strength.  Yet  in  other  cases  grooves  and  flat  spots  would  considerably 
increase  the  danger  of  buckling.  For  a  refined  treatment  of  imperfection 
effects  it  would  be  necessary  to  know  the  effect  of  a  variety  of  sizes, 
shapes,  and  positions  of  blemishes.  Also  it  would  be  necessary  to  know 
the  mode  shapes,  both  symmetric  and  unsymmetric.  Finally,  it  is 
postulated  that  the  position  of  the  imperfection  with  respect  to  troughs 
in  the  deflected  shape  would  have  an  effect  on  buckling  pressures.  Such 
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complex  relationships  make  some  sort  of  minimum  load  criterion  very 
welcome  although  somewhat  Improbable. 

Testing  offered  a  further  possibility  for  imperfect  conditions. 

The  latter  can  change  from  test  to  test  while  rim  slope  and  surface 
imperfections  will  not.  Hence  variations  in  buckling  pressures  for 
a  single  shell  are  probably  due  to  changes  in  test  conditions.  One 
source  of  error  is  the  accuracy  of  clamping  the  shell  into  the  test 
chambers.  If  the  shell  were  placed  slightly  off  center,  the.  shell 
would  not  be  held  correctly  and  symmetrical Ij^  but  in  some  unsymmetrical 
manner.  This  would  tend  to  Induce  unsymmetric  buckling  and  consequently 
lower  loads.  The  variation  in  buckling  pressure  which  may  be  attributed 
to  this  cause  is  about  7^  as  can  be  computed  from  the  data  of  Table  1. 
During  the  air  pressure  tests,  the  nature  of  the  effect  of  clamping 
conditions  was  studied  for  one  shell.  It  was  found  that  the  buckle  could 
be  made  to  form  in  any  part  of  the  shell  by  shifting  the  shell  around 
slightly  in  the  test  chambers.  Early  in  the  testing  program  this 
importance  of  clamping  conditions  was  recognized  so  that  great  care 
was  exercised  in  this  regard.  In  the  preliminary  tests  when  no 
special  precautions  were  taken,  the  reduction  due  to  Incorrect  edge 
conditions  was  as  much  as  50^. 
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-  Comparison  with  Other  Work  - 

The  buckling  loads  of  this  report  seemed  to  fall  in  a  sort  of 
middle  range  between  those  of  Tsien  and  Kloppel  and  Jungbluth,  but 
are  quite  near  the  results  of  Kaplan  and  Fung,  and  Homewood,  Brine, 
and  Johnson.  This  was  some  verification  that  the  results  from  the 
plastic  material  were  valid. 

In  the  History  it  was  pointed  out  that  if  the  energy  criterion 
of  Tsien  were  correct,  buckling  loads  would  depend  on  the  ambient 
energy  level  and  not  be  very  reproducible  from  test  to  test.  However, 
as  was  mentioned,  the  values  of  are  very  repeatable  under  unchanged 
conditions.  The  criterion  was  also  used  to  predict  that  buckling  loads 
under  constant  pressure  would  be  about  half  those  under  constant  volume 
testing  conditions,  Kaplan  and  Fung  made  both  types  of  tests  and  found 
no  indication  of  a  difference;  however,  their  experiments  were  in  shallow 
shells  beyond  the  range  of  Tsien 's  criterion.  The  tests  of  this  report 
are  evidence  that  the  constant  pressure  buckling  load  may  actually  be 
above  the  constant  volume  load.  Such  a  circumstance  is  not  explainable 
either  by  the  energy  criterion  or  by  the  usual  maximum-point-on- the- 
cu  ’'e  criterion. 

Tlie  main  importance  of  this  experimental  work  was  in  the  repeated 
testing  and  in  the  study  of  effects  which  could  be  investigated  by 
this  type  of  testing.  Since  the  other  investigators  used  metal  shells 
which  usually  yielded  on  buckling,  only  one  critical  load  was  obtainable 
and  many  questions  were  raised  as  to  why  the  shell  buckled  at  that 
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prebsure,  Soms  of  the  possible  causes  have  been  studiefi  in  this  report 


by  using  shells  which  buckled  elastically  and  could  be  re-buckled  any 
number  of  times.  Thus  the  importance  of  clamping  conditions  has  been 
pointed  out  and  the  time-dependence  of  buckling  loads  in  material  which 
creeps  was  considered.  It  was  possible  to  conclude  that  the  buckling 
was  due  to  a  geometric  instability  and  not  to  material  yield  or  to  some 
local  failure.  Since  the  shells  were  buckled  many  times  it  was  discovered 
that  the  dimple  could  appear  in  many  positions  in  the  shell  and  that  the 
buckling  loads  and  buckle  positions  were  related.  Finally,  an  unbuckling 
curve  was  determined  which  may  shed  further  light  on  the  buckling 
process . 
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CHAPTER  5  ANALYTICAL  STUDY 


A  Introduction 

The  main  purpose  of  an  analysis  in  a  report  on  buckling  is  the 
determination  of  buckling  loads.  However,  with  the  large  disparity 
between  analytical  results  of  previous  investigators  and  test  results, 
there  was  not  much  hope  for  finding  better  results  for  from  this 
analysis.  Rather,  the  computations  were  made  to  study  the  nature  of 
the  load-displacement  curves.  The  loading  curves  were  plotted  both 
as  load  versus  center  displacement  and  load  versus  volume  change. 

In  addition,  it  was  desired  to  know  the  variation  of  energy  along 
these  curves  in  crder  to  assist  in  a  consideration  of  Tsien's  energy 
criterion , 

The  technique  of  calculation  was  the  energy  or  variational 
method.  Since  this  approach  provides  a  series  of  equilibrium  posi¬ 
tions,  the  necessary  curves  could  be  plotted  from  the  coordinates 
(pressure  and  deflection)  of  these  positions.  The  required  equation 
of  total  potential  energy  was  first  derived  from  the  familiar 

principle  of  virtual  work.  Nonlinear  strain-displacement  relations 

(22 ) 

were  taken  from  Love  for  use  in  the  energy  equation.  Finally  the 
energy  equation  was  simplified  to  contain  only  terms  important  for 
large  axisymmetric  deflections  of  deep  spherical  shells  with  a 
buckle  at  the  apex.  Most  previous  derivations  either  were  restricted 
to  shallow  shells,  used  different  strain-displacement  relations,  or 
were  linear  and  not  applicable  to  the  large  deflection  problem. 
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In  ora<_r  to  continue  the  an«ilysis,  some  estimate  had  to  be  made 


of  an  equation  for  deflected  shape.  The  usual  procedure  at  this 
point  in  the  problem  has  been  to  propose  a  series  of  functions  whicu 
seemed  to  fit  the  equations  mathematically.  The  estimate  used  in 
this  investigation  was  taken  from  observation  of  the  actual  buckling 
tests.  The  deflected  shape  was  assumed  to  be  a  sum  of  tenns :  one  for 
the  uniform  deflection  as  predicted  by  small  deflection  theory,  one 
for  the  apparent  shape  of  the  buckle,  and  several  others  to  satisfy 
boundary  conditions.  The  formula  represents  quite  well  the  appearance 
of  the  deflected  shape,  particularly  in  the  post-buckling  range. 

The  first  maximum  of  the  load-del lection  curve  was  assumed  to 
be  the  buckling  load.  This  maximum  point  was  found  graphically  by 
plotting  a  number  of  cquilibr iu.ii  positions  near  the  maximum. 

The  calculations  were  carried  out  partially  on  the  IBM  7090 
computer  because  of  the  complexity  of  terms  in  the  energy  equation. 
Equilibrium  positions  were  found  foi  unrestricted  deflection  and  for 
deflection  where  the  buckle  increased  in  depth  but  the  average  de¬ 
flection  over  the  shell  remained  constant. 
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B  Derivation  of  Equations 


In  this  section  the  well-known  principle  of  virtual  work  was  used 
to  derive  the  principle  of  minimum  total  potential  energy.  In  the 
application  of  the  latter,  a  functional  known  as  the  total  potential 
energy  was  calculated  for  use  in  determining  equilibrium  equations. 

In  deriving  the  equation  for  the  functional,  various  "order  of  im¬ 
portance"  studies  were  made  on  the  terms  in  the  equation  so  that 
negligible  terms  could  be  discarded. 

The  principle  of  virtual  work  may  be  stated  as  follows  (see 
K.  Washizu,  ref.  21 j ;  If  infinitesimal  virtual  displacements  of  a 
body,  u,  V,  w,  which  are  compatible  with  the  prescribed  boundary 
conditions  are  imposed,  then  at  an  equilibrium  position,  the  principle 
of  virtual  work  gives 


Volumt 


"//(xSu  y  ZS<Ar)c/S  =  O  (s-i) 

Sur  face. 

The  quantioies  X,  Y,  Z  are  forces  on  the  surface  acting  in  the 
coordinate  directions,  and  u,  v,  and  w,  are  displacements  in  the 
same  directions.  The  strains ,  ,  <£22  t  ^ as  >  ^ 2 3  >  3  ■  &nd  £^2 

are  derived  from  the  displacements.  Also  the  stresses  are  specified 
as  functions  of  the  strains  by  stress-strain  relations,  such  as 
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From  the  principle  of  virtual  work  Washizu  has  derived  the 
principle  of  minimum  potential  energy  which  may  be  stated  as 
follows:  Of  all  the  admissible  u,  v,  and  w  which  satisfy  the 
prescribed  geometric  boundary  conditions,  the  actual  displacements 
are  given  by  the  extremum  condition  of  a  functional^  defined  as 


\/oluinne 


t  9'ir  -f  Z  ur)  0I S 

Surfuce 


(5-2) 


If  the  first  variation  of  TT.  the  total  potential  energy,  is  taken 
with  respect  to  the  displacements  and  strains,  Eq .  (5-1)  is  obtained. 
The  same  minimization  procedure  is  valid  for  any  number  of  independent 
displacement  parameters  and  is  not  restricted  to  three  orthogonal 
components . 

In  the  usual  application  of  the  energy  method,  the  total 
potential  energy  is  evaluated  and  then  varied  with  respect  to  the 
displacement  parameters.  The  first  vai ration  provides  the  minimum 
condition  ^ 

STf  -  t  SlT  t  f  ”  ’  -  O 


ov  r  T — 
Jtr  Jur 


Since  the  displacements  6u,  Sv,  5w,  .  .  are  independent  and 
may  be  specified  arbitrarily,  each  partial  derivative  must  be  zero 
to  ensure  that  Sir-  0.  Therefore  the  minimization  and  equilibrium 
conditions  are 


••• 


Ju. 


Jv 


(5-3) 
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If  the  displacement  components  are  not  independent,  a  slightly 
modified  procedure  can  be  followed.  For  a  two-variable  problem  the 
usual  equilibrium  condition  is 


Srt  =  ^Su  7- 

du  Jyr 


(5-4) 


but  if  u  and  v  are  not  independent,  there  is  a  constraint  condition 
in  the  form  M(u,v)  =  0.  Therefore  it  is  also  true  that 


SAI  =  =o 

dU  dV 


or 


Replace  this  expression  into  (5-4)  and  obtain 

J(Z 

Since  Sv  can  now  be  varied  arbitrarily,  the  equilibrium  condition  is 

» 

found  to  be 

Jv-Ju  Ju.  "Sv- 

For  convenience  the  total  potential  energy  is  often  separated 
into  two  parts ,  thus 

77^  =  U  -  W,  (5-5) 


where  U  is  the  strain  energy  and  W  is  the  work  of  the  external  forces. 

The  present  problem  may  be  considered  as  a  linearly  elastic 
case  of  plane  stress,  that  is,  the  stress  (T'j,  is  taken  as  zero.  Also 
Kirchhoff  s  assumption  that  a  line  in  the  shell  which  is  normal  to 
the  middle  surface  before  deformation  remains  normal  after  deformation, 
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will  be  taken.  This  merely  requires  that  the  shear  strain  will  be 


neglected  in  finding  deflections.  Further  it  is  assumed  that  it  is 
reasonable  to  neglect  terms  which  appear  to  be  small  accord.^ng  to 
some  consistently  applied  criterion.  In  this  problem  the  strain 
energy  takes  the  simple  form 


Volume. 


and  the  stress-strain  relations  are 

^*'4  ^  ) 

^ 5~7  ) 

The  quantities  e^n,  *^22  dimensionles.s  curvatures  of  the  middle 
surface  in  the  0  and  0  directions  and£ji,  C  22  strains  of  the 

middle  surface  in  those  two  directions. 

Substitution  of  the  relations  (5-7)  into  Eq.  (5-6)  and 
integration  in  the  direction  of  0  gives 

(5-8) 

The  required  strains  and  curvatures  are  taken  from  the  nonlinear 

strain-displacement  equations  of  A.  E.  H.  Love'"^’.  The  following 

(23a) 

procedure  is  outlined  in  Novozhilov  .  The  relations  are  specialized 
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here  for  the  case  of  symmetrical  deformation  of  spherical  shells. 


The  strains  are 

^  ~^zf] 

<?3  *  "^^3  ^ 

=  ^3  ^  >)  '  <^33  ^ 


where 


/ 

a¥  z 


C3 


/ 


<fy2 


(u^cof^  T>-  CJi) 


L—  ( ^  )  7^  f^C 


^53  == 

=  d^'c _ L _ (  dfic  -  ) 


(5-9) 


(5-9) 


Here  u^  and  w^  are  displacements  tangential  and  normal  to  the 
middle  surface  and  z  is  the  distance  outward  from  the  middle 
surface. 

In  order  to  proceed  further  it  is  necessary  to  find  a  criterion 

to  determine  the  relative  importance  of  displacement  terms  so  that 

small  terms  may  be  discarded.  Edge  loads  and  point  loads  on  shells 

both  produce  sudden  changes  in  geometry  and  high  moments  just  as  the 

(24 ) 

buckle  does.  From  the  analysis  of  point  loadings  by  Reissner'  and 

(23b) 

edge  loadings  by  Novozhilov  the  relative  importance  of  the 

deflection  terms  can  established.  In  these  analyses,  u^ ,  the  tangential 
deflect!  '’omponent,  is  small  with  respect  to  w^ ,  the  normal  component. 
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In  other  words,  w  is  of  the  order  of  u  times  m,  where  m  is  a  number 
c  c 

which  is  larger  than,  say,  10.  Also  differentiation  by  <f)  or  multiplica¬ 
tion  by  cotjzi  raise  the  order  by  one  step,  that  is,  increase  the  magnitude 
by  one  power  of  m.  Therefore  u^cot0  and  u^  are  of  vne  order  of  (u^)(m), 
w'  and  u"  are  of  the  orderof  (u^)(m^),  and  w^'  and  w^cot0  are  of  the 

order  of  (u  ) (m*  ) .  These  orders  of  magnitude  apply  specifically  near 
c 

the  apex  where  cot0  is  m".ch  greater  than  one  and  are  no  t  as  applicable 
at  the  edge  of  a  deep  shell.  In  this  analysis  it  is  assumed  that  the 
highest  order  nonlinear  terms  in  the  strain-displacement  relations  are 
comparable  to  the  highest  order  linear  terms  and  that  all  other  linear 
and  nonlinear  terms  are  negligible. 

From  Kirchhoff's  assumption  it  follows  that  ^13  must  vanish  and 
the  plane  stress  assumption  requires  that  =  0.  Therefore 


(5-12) 


Le 


et  u  =  u  + 


w^  =  w  ^ 

where  u  and  w  are  displacements  of  the  middle  surface  and 

and  X  33"©  linear  variations  of  displacement  through  the  thickness. 

Substitute  these  values  into  the  equations  for  £^3  and  £33  and 


(5-13) 
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the  results  are 


'  =  ( I -tXf  i- ■0’' 


Let  e^  =  u'  +  w 
e^  =  ucot0  +  w 


e,  =  w'  -  u 


(5-14) 


where  u  and  w  are  dimensionless  terms  for  the  displacements  of 
the  middle  surface,  that  is.  u  ,  w  =  -f  .  The  primes  refer  to 
differentiation  with  respect  to  0,  Then 

^  (5-15) 

X=  -J-e/ 

Now  return  to  £  and  Use  the  substitution  (5-13)  and  let 

4  ^  ^  r ‘^/z  ^  Then 

■=  ~  ur  '  ■& uuw' t y(Ari/>r"  ~  (5-16) 

-  ^  coi-<^  r*  ^  ~  uur'  cof^  ^ 

=  -CATcoff^ 

The  results  above  have  bei-n  found  by  keeping  only  the  highest  order 

linear  and  nonlinear  terms  as  determined  by  the  order  of  magnitude 
analysis . 
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The  solution  of  the  problem  is  greatly  facilitated  by  solving 
for  u  in  terms  of  w  so  that  only  w  remains  as  an  independent  variable. 
The  necessary  condition  is  found  by  minimizing  the  total  energy  with 
respect  to  u.  Since  u  appears  only  in  the  extensional  energy, 


STT  -  SU  ^  ^ 


which  is  an  equilibrium  condition  for  a  homogeneous  case.  This 
may  also  be  transformed  to 

=  -Oiu)uT'- {rr'iAr” -  cof^ 


The  fourth  tenn  is  of  smaller  order  than  the  others  in  the 
equation  and  may  be  omitted.  The  rest  can  be  arranged  for  in¬ 
tegration  as  follows: 


~  iLz^^ cof/ 


(5-18 ) 
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[ntegration  once  gives 


y- i/cdf 9^  ~ 


a.  /*  2 

—  -(/  ■yu)(4r  -  “*  I  -y  ^ 


(5-19) 


Here  it  is  well  to  notice  that  there  is  a  simple  expression  for  the 
siuu  of  the  extensional  strain 


4  5,:  = 


(5-20) 


Notice  that  if  v  is  a  unifora  expansion  this  equation  is  incorrect 


since  i//  ^2  ~  for  Cat  ar  co^stonf. 

Another  integration  of  (5-18)  gives 


(5-21 ) 


'  Cofcf<9/(6  ^9^  y  /fO-cos^)  ^  3 

7  r  -X— — 


Imposition  of  the  boundary  conditions  u  =  0  at  0  =  0,o( 
serves  to  specify  A  and  B. 
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Since  the  sum  of  extensional  strain,  £  n  +  ^*2.  (see  5-20) 
is  easy  to  calculate  but  the  product,  £>  ^  strain 

energy  formula  will  be  modified,  using  Bq.  (5-7).  The  extensional 
energy  is 


The  last  integral  can  be  integrated  partially  using  equation 
(5-17). 


The  first  term  is  zero  at  the  lower  limit  and  only  u'  is  nonzero 
at  the  upper  limit.  Therefore  the  integral  of  the  strain  product 
becomes : 

-  £'J7  7r'a^s/»o(  7^/70 

•f’ 
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The  integral  is  of  the  order  of  l/m*  times  the  evaluated  term  and  can 
therefore  be  integrated  approximately  or  neglected. 

The  forinule  for  bending  energy  can  also  be  transformed  by  partial 
integration  to  a  simpler  form. 


The  second  term  in  is  integrable  since 

D 

idf^cospij  y  S//f  pi  c/d 

The  evaluated  term  is  identically  zero  and  the  integral  is  of  the 

order  of  l/m*  with  respect  to  the  first  term  in  U.  .  The  nonlinear 

b 

terms  in  the  bending  energy  are  found  to  be  second  order  with  respect 
to  the  comparable  extensional  energy  terms  and  are  therefore  omitted. 

The  work  done  can  be  considered  as  a  scalar  product  of  force  and 
distance,  thus 


W  ^  fjj( a  COS )(  -t  ((T  cos2 ) 


(5-22) 


where  cos  X.and  cos  Z  are  cosines  of  the  angles  between  the  surface 
normal  and  the  0  and  Z  coordinate  directions.  Here 


CoS  X  =  ^ 

cosZ  -  !  -f-  X 


(5-23) 
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Therefore 


W -  -  ue^ S/'r? ^ 


or  approximately 


IV  -  P7na^<^jurS/h^ 


(5-24 ) 


The  energy  equation  can  now  be  written  in  the  form: 


ft  ^  C(,i^  t  -dL 


whore  the  dimensionless  quantities  are 


ff  =  TTC/-^') 

£^h7ra^ 

u,  =  itj±  f 'l,  ^  V/7/  y/ 


/ 


S//7o(  7^/7<?<  6/ 


^6 


__  / 


/  -✓ 


(s.. 

Cos^pi 


-  .  ^  ..  S'/V? ^ 


Q 

L 


oi 

IAT  sif^^  ci ^ 


(5-26  ) 
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C  Deflected  Shape 


The  deflected  shape  to  be  introduced  into  the  equation  for  the 
energy  functional  is  probably  the  main  difference  between  this  analysis 
and  previous  work.  In  the  past  the  deflection  was  taken  in  a  series 
w’ th  as  many  as  200  unknown  coefficients.  With  fewer  terms  the  shape 
was  too  restrcted  and  poor  results  were  obtained.  Also  the  number 
of  terms  increased  rapidly  as  increased.  For  a  deep  shell  (large  ^  ) 
a  power  series  might  prove  hopeless. 

In  this  analysis  a  deflected  shape  was  chosen  which  seemed  to  fit 
the  appearance  of  the  shell  surface  during  loading.  Up  to  buckling 
the  deflection  seemed  to  be  uniform  over  most  of  the  shell  (this  was 
observed  under  air  pressure  where  loading  could  be  applied  rapidly  so 
that  deflection  was  visible).  At  buckling  a  dimple  formed  which 
was  a  shallow,  bowl-like,  circular  depression.  The  edges  seemed  to 
be  fairly  well  defined  due  to  the  reflection  of  light  from  the 
region  of  sharp  curvature  at  the  rim  of  the  depression.  Therefore  it 
was  possible  to  speak  of  a  diameter  or  size  of  the  buckle.  This 
sharpness  of  curvature  indicates  the  presence  of  large  and  rapid 
variations  in  shape,  shear,  slope  and  moment  in  limited  areas  of  the 
surface.  Such  regions  of  high  derivatives  are  difficult  to  depict 
analytically  and  require  a  long  series  of  terms  to  make  an  adequate 
fit.  In  the.- she  11  -the  buckle  appeared  as  almost  a  separate  region 
in  an  otherwise  unchanged  spherical  surface.  No  wonder  many  of  the 
first  investigators  analyzed  the  buckle  region  as  a  separate  entity. 
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The  buckles  seen  in  the  tests  were  symmetrical  only  for  the  shallow 
shells  with  one  exception.  However  the  differences  in  buckle  position 
apparently  caused  fairly  minor  changes  in  Therefore  if  buckling 

loads  for  axisymmetric  buckles  were  determined,  the  load  associated  with 
unsymmetric  buckles  should  be  just  5^  or  10^  lower.  The  axisymmetric 
case  only  is  considered  here  because  it  is  simpler  mathematically  and 
should  depict  adequately  the  nature  of  the  loading  curve. 

Instead  of  trying  to  depict  the  whole  shape  of  the  shell  surface 
as  one  term  or  one  series  of  terms,  the  shape  was  divided  naturally 
Into  its  two  components:  the  uniform  deflection  and  the  buckle  de¬ 
pression.  No  mode  shapes  were  considered  since  none  were  seen  on  the 
shells.  These  small  undulations  in  the  surface  were  probably  just  too 
small  to  be  detected  visually  without  special  equipment.  Omitting 
modal  deformations  from  the  assumed  deflected  shape  probably  led  to 
a  somewhat  higher  buckling  load. 

The  equation  for  deflection  which  was  chosen  was  as  follows: 

-  c/r^[str)c(oi-^)  /  cos 
-ur,  (s//t  /  cos  6^ )c~ 

(5-27) 

^ot  /  CoS  Lot.  J(^~  ^ 

Ceos  fi  -cos-ot) 

S//>ot 

Each  term  here  had  a  specific  purpose.  The  first  term,  Wj,,  is  a 
uniform  expansion  and  v in c(o( -  i-cOS  cO>t^p^  ^  ^ovides  for  the 
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claraped  edge  condition.  Both  of  these  deflection  terras  are  known  from 
small  or  linear  rtheory.  The  third  terui,  Wj(.sin  b0  +  cos  h(^)e 

describes  the  shape  of  the  dimple  which  appears  at  buckling.  For  con¬ 
venience  these  first  three  terms  were  plotted  in  Fig.  4B  in  various  ways 
Notice  that  the  third  term  has  no  slope  at  0  =  0,  and  that  it  goes  to 
zero  at  b0  =  2.35  radians.  For  larger  values  of  b0,  the  contribution 
of  this  term  becomes  very  small,  l/b  is  a  measure  of  buckle  extent 
since  the  diameter  of  the  dimple  is  about  4.7a/b. 

The  succeeding  two  terms  were  needed  to  satisfy  boundary  condi'^ions 
for  w  and  w'.  With  the  clamped  edge  it  was  necessary  that 

w~w'=0  at0=:C)(,  the  edge 

w'=  0  at  0  =  0,  the  apex 

and  w  had  to  be  finite  at  0  =  0.  These  boundary  conditions  are  not 
satisfied  exactly  but  approximately,  since  the  edge  correction  has 
both  a  deflection  and  slope  at  the  apex.  The  inflection  in  the  cen¬ 
ter  of  the  buckle  is 

'*'o  ~  fsin  bo<  +  cos  b</  )e 

-bo< 

bWjSin  bo^  e  sino^ 

but  for  most  purposes  it  was  adequate  to  assume  that  the  buckle  depth 
at  the  apex  was 

-  wo  +  Wj 

(the  terms  were  written  with  different  signs,  but  in  actuality  both 
are  inward  deflections  and  therefore  edd  since  Wj,  is  always  negative.) 
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In  general  both  c  and  b  will  be  large,  that  is,  between  10  and 
100.  This  means  that  the  number  of  important  terms  in  the  series 
expansion  for  w  would  be  extremely  large  and  this  seems  to  be  the 
case  for  deeper  shells.  The  factor  b  is  the  basis  for  all  order  of 
magnitude  analyses  and  exactly  fits  the  scheme  earlier  specified. 

Near  =  0,  it  is  assumed  that  b0  is  of  the  order  of  unity.  This 
deflection  equation  is  reasonable  for  shells  of  a  certain  depth 
only,  where  the  edge  and  buckle  terms  do  not  significantly  overlap 
and  is  not  necessarily  applicable  to  very  shallow  shells.  For  the 

— C  oC 

deeper  shells  it  is  possible  to  let  e  equal  zero  in  the  evaluations 

-bc< 

of  theintegrals .  e  cannot  be  neglected  since  it  increases  on 
differentiation  with  respect  to  b. 

The  energy  equation  can  be  written  in  the  following  general 

form: 

/  ^  ^03  007^ 

(5-28) 

/  ^  -  ^0/ 


The  E's  contain  exponential  and  trigonometric  functions  of  b  and  c. 
Equilibrium  equations  are  now  found  by  calculating 

Si>  de¬ 


ls 


The  problem  is  simplified  by  using 


(^Q  * 


/- 


C  = 


(5-29) 


which  can  be  obtained  by  evaluating  0  and  0,  neglecting 

au>1  ^ 

terms  in  Wj .  This  merely  means  that  w-  is  assumed  to  be  the  correct 


linear  deflection.  But  if  this  is  so  then  the  terms  with  coefficients 


Eji  and  must  also  cancel,  as  Friedrichs  has  pointed  out,  and  could 
be  omitted  from  the  energy  and  equilibrium  equations.  That  this: was. 


the  correct  procedure  can  be  demonstrated  using  a  simple  potential 
energy  formula  which  contains  both  a  linear  and  nonlinear  term. 

Let  TTo  —  f^S  -  S  ^2  (5-30) 

^  can  be  considered  linear  since  its  derivative  is  linear.  This 
special  definition  applies  since  TTis  actually  the  integral  of  the 
physical  quantities  which  are  linear  or  nonlinear.  Also  the  deflection, 
$  ,  is  broken  into  a  linear  term  which  is  valid  near  5=0,  and  a  non¬ 
linear  term,  S,  . 

S  -  So  S, 


Pi,  -  C-  -C  iX-i  C 
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Notice  that  if  Wq  =  So  >  *1  ~  /Ti  then  Eq.(5-31)  is  very 

similar  to  Eq.  (5-28).  The  former  is  much  simpler,  however,  and  can 
be  studied  to  determine  the  importance  of  the  terms.  The  general 


solution  can  be  found  from 

In  this  case  there  is  a  maximum  or  buckling  pressure  which  occurs 
at  df*  =■  >6/  -  ^  £3  £  -  O 
or  £  -  and  P  =  - 

Now  the  same  result  could  be  found  using  S~Se'^^/»  When  Eq.  (5-31) 


is  differentiated  with  respect  to  St,  the  following  is  obtained: 


-  -- -ci  -ax  -  p 


(5-32) 


and  P  occurs  at  s  —  ^ 

max 


(5-33) 


zi/ 


which  is  the  result  obtained  before  for  the  deflection  at  buckling. 
Notice  that  in  7^ the  terms  which  contribute  to  (5-33)  are  only  ^ 

-3)fa  S  o%\,  and  -^2  S  1  ’  Since  So  ts  the  linear  solution, 

$0  =  P/ki .  Using  this  result  and  the  deflection  at  maximum  P  in 


Eq.  (5-32)  gives 


^  _  // 


In  finding  P  ,  four  terms  in  /TV  are  needed.  Besides  the 
max 

three  which  serve  to  determine  the  critical  deflection,  -3kj  f 


is  required.  If  the  latter  is  omitted,  the  value  of  P  is  twice 

max 
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the  correct  value.  The  interesting  result  is  therefore  that  the 
linear  solution,  So>  roust  be  included  in  two  nonlinear  terms,  8^ 
and  s.r  in  order  to  give  the  correct  results.  In  passing  it  should 
be  noted  that  the  terms  in  TT^  containing  Eq,  So^i»  ^Sqi 
pSi  cancelled  out  since  So  the  linear  solution. 

When  the  results  of  the  preceding  study  were  applied  to  the  terms 
in  fT,  Eq.  (5-28),  it  was  found  that  the  terms  with  coefficients 
Eao>  Sill  Sjo,  and  could  be  dropped  but  that  all  the  others  had 
to  be  kept.  Finally  the  energy  equation  took  the  form 

"TT"  “■  I  ^ m  "f"  *”  ^oi  ^  ^oy  ^  (5—34  ) 


(5-35) 
(5-36) 
(5-37  ) 


where  the  D's  derivatives  with  respect  to  b  of  the  E's.  The 
E's  found  in  Eq .  (5-34)  are  given  in  Appendix  C. 
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D  Solution 


The  procedure  for  solution  of  the  problem  is  to  satisfy  the  two 
equations  (5-36)  and  (5-37)  simultaneously  for  Wj  and  b  at  various 
values  of  Q,  the  applied  pressure  parsuneter.  Each  set  of  values  for 
Q,  Wj  ,  Euid  b  which  satisfy  the  equations  define  an  equilibriiun  posi¬ 
tion  for  the  shell.  If  a  series  of  quantities  for  Q  and  Wj  from 
such  positions  are  plotted,  a  loading  curve  is  obtained.  In  this 
report  it  was  found  convenient  to  plot  -w^,  versus  -w^  +  Wj  since 
-Wfl  is  proportional  to  pressure  and  -Wj,  +  w^  is  the  central  deflection. 

The  first  attempt  at  a  solution  was  made  with  a  somewhat  less 
TY  — b0( 

complete  formula  for  // .  Since  e  is  rather  small  for  large  b, 
all  terms  with  an  exponential  were  excluded.  This  meant  that  all 
trigonometric  functions  were  also  omitted  since  they  occur  only  in 
conjunction  with  the  exponential.  Hence  the  E's  contained  only 
powers  of  b,  making  the  equations  rather  easy  to  solve.  Terms  con¬ 
taining  e  ^  were  also  neglected  in  the  expression  for  w,  Eq.  (5-27). 
And  finally  E2iWoWj  was  neglected  because  it  did  not  seem  reasonable 
for  the  linear  solution,  w^,  to  be  of  imp>ortance  in  nonlinear  terms. 

The  two  equilibrium  equations,  0  and  0,  were  solved  for  a 

number  of  values  of  Q.  Eauilibrium  positions  were  obtained  first  far 
into  the  post-buckling  range,  and  the  solution  was  then  continuevi, 
working  backwards  towards  buckling.  The  post-buckling  part  of  the 
curve  was  about  right.  However,  at  a  point  near  what  should  have 
been  the  maximum  point  of  the  curve,  the  pressure  rose  sharply 
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towards  infinity  and  the  value  of  b  became  imaginary.  Although 
this  solution  atteir.pt  failed,  it  gave  some  useful  information 
for  assistance  on  the  next  try.  It  was  found  that  b  would  have  a 
value  near  that  of  c  and  that  the  value  of  b  rose  to  a  maximum 
near  P  .  Thus  a  curve  of  b  versus  deflection  should  be  similar  to 


cr 


a  load-deflection  curve, 


The  difficulties  experienced  in  the  first  try  could  be  traced 

-bey  -hoi 

mainly  to  the  omission  of  all  terns  containing  e  .  True,  e 
was  much  smaller  than  l/b,  with  which  it  was  usually  compared, 
especially  if  b  were  approximately  equal  to  c.  However,  when  the 
derivatives  of  both  terms  were  taken,  the  situation  was  different. 
Then  the  comparison  was  between  be  and  l/b^  and  the  two  terms  are 
nearly  equal  for  b  =  36 . 

In  the  second  attempt  at  a  solution,  the  procedure  outlined 
in  most  of  this  chapter  was  used.  The  nonlinear  terms  Egi^^Wj  and 


'^■^re  retained,  all  quantities  containing  e  ,  e  ,  e 
and  e  were  retained,  and  £20'^!  >  0 ‘*'o  were 

cancelled  out.  The  remaining  E’s  were  very  complicated  functions 
of  b  as  can  be  seen  from  the  formulas  in  Appendix  C.  The  normal 
procedure  for  the  solution  of  Eq.  (5-36,  5-37)  is  to  specify  Q 
and  solve  for  w^  and  b.  However  it  is  also  possible  to  specify 
b  and  solve  for  Q  and  Wj  if  that  is  more  convenient. 

In  the  present  case  computations  were  made  by  evaluating  the 
E's  and  D’s  at  certain  values  of  b.  b  was  successively  assigned 


-be*  -2bo(  -3bcx 
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values  for  each  percent  of  c,  that  is.  b(i)  =  O.Ol(i)c  where  i  was 
an  Index  running  from  30  to  120.  This  range  of  b  was  suggested  by 
the  results  of  the  first  attempt  at  solution  and  the  interval  was 
chosen  to  give  reasonable  results  for  the  derivatives.  The  latter 
were  computed  by  a  four  point  central  difference  formula  from 
reference  (25).  For  example 


and  here  the  range  of  i  was  narrowed  to  32  to  118 . 

The  computations  of  the  E’s  and  D's  were  programmed  in  Fortran 
language  for  use  on  the  IBM  7090  at  M.I.T.  The  program  was  written 
in  essentially  the  order  given  for  the  calculations  in  Appendix  C. 

When  these  evaluations  had  been  made,  the  simultaneous  equations 
of  equilibrium,  (5-36,  5-37)  were  ready  for  solution.  Since  each  set 
of  these  equations  were  for  a  fixed  value  of  b,  Wj,  and  Wj  (or  Q  and  Wj  ) 
were  treated  as  unknowns.  The  remaining  calculations  were  made  by 
hand.  The  method  used  v.-is  to  try  to  solve  a  few  sets  of  equations 
for  b  near  0.9  c.  When  a  solution  of  any  set  was  obtained  then  a 
set  for  a  nearby  value  of  b  was  attempted.  To  aid  in  solution,  the 
two  equilibrium  equations  were  rewritten  as 


- 


(5-38) 
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Then  the  equations  were  solved  by  relaxation.  That  is,  values  for 
Wq  and  were  estimated  and  placed  in  (5-38)  and  a  new  value  for 
Wq  was  found.  Then  with  the  estimate  for  Wj  and  the  new  value  for 
Wg,  Wj  was  evaluated  using  the  second  equation.  The  process  was 
repeated  until  a  set  of  solutions  was  found  or  until  it  became 
apparent  that  the  values  of  one  or  both  of  the  unknowns  would  lie 
outside  the  range  of  interest.  When  successive  values  of  Wq  and 
agreed  to  three  figure.**,  the  equations  were  considered  solved.  The 
values  of  and  Wj  obtained  were  plotted  immediately  to  assisv  in 
estimating  the  next  solutions. 

For  the  graphs  shown  in  Fig.  16,  the  plot  was  made  of  -Wg  versus 
-Wg  +  Wj  and  then  the  ordinate  scale  was  changed  to  make  a  pressure 
versus  center  deflection  curve.  The  pressure  versus  volume  change 
curves  were  made  by  plotting  -Wg  versus  L  as  given  by  Eq .  (5-26). 

The  preceding  calculations  were  all  for  "free"  deflections .  that 
is  Wi  and  b  were  varied  independently  to  minimize  the  total  poten¬ 
tial  energy.  However,  when  an  experiment  was  run  with  volume  control, 
this  freedom  was  not  present.  To  duplicate  this  case  analytically, 
a  restriction  was  applied  to  w^  and  b  so  that  the  volume  L  could  be 
held  constant  while  deflection  and  pressure  changed.  The  method 
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used  here  was  outlined  in  th«  discussion  following  Eq .  (5-4). 

M(u,v)  =  0,  the  restraint  condition,,  was  replaced  by  L  -  L®  = 

where  !•  the  value  of  L  at  some  designated  point.  For  instance, 

the  value  of  L  at  P  was  set  equal  to  L^.  Then  with  th«  constraint 

cr 

condition,  equilibrium  positions  shown  in  Fig.  15  duplicated  the 
curve  of  an  actual  shell  undergoing  constant  voluHe  testing  (see 
Fig.  10). 

For  both  free  and  volume  control  testing  It  was  possible  to  es¬ 
tablish  cux^es  of  energy  vereus  deflection.  These  were  simply  made 
by  uslnc  the  values  w^ ,  b,  and  Q  which  were  already  known  at  specific 
points  and  evaluating^  from  Eq.  (5-28). 

During  the  solution  of  the  equatle&s,  It  was  found  that  the  terms 
with  Bi2  **^<1  ®3i  were  very  important  in  the  initial  portion  of  the 
loading  curve  up  to  buckling.  Also  in  this  Initisl  portion  the  terms 
with  £o3  and  E04  were  not  required  end  Eq,  wms  contrllMtlng  s  very 
small  amount  even  st  buckling.  The  situation  was  rersrsed  for  equili- 
briom  positions  in  the  large  deflection  range.  The  terms  with  Eq, 
and  Eo4  were  dominant,  and  Bj,  and  Ejx  were  essentially  negligible. 

The  linear  term  Bqb  Importance  throughout  the  range  of  para¬ 

meters  . 

One  effect  of  the  unimportance  of  Ej^WqW^  and  IjbWqW^  for  large 
deflections  was  the  poor  definition  of  w.^  and,  therefore,  pressure. 
Thus  when  w^  was  large,  one  equation  was  solved  for  almost  in¬ 
dependently  of  Wq.  Thus  the  usual  Interaction  of  simultaneous  equa¬ 
tions  was  not  present  aod  the  weiues  of  pressure  were  somewhat 
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doubtful.  Also  for  large  deflections  the  potential  energy  was  not  well 


defined  since  it  was  a  small  difference  between  large  quantities.  Since 
the  term  was  the  largest  energy  term  as  w^  became  greater  than, 

say,  0.01,  it  seemed  apparent  that  terms  in  wf,  wf ,  etc.  would  be  of 
importance . 
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E  Analytical  Results 


A  detailed  analysis  as  outlined  in  the  Solution  Section  was 
completed  only  for  a  configuration  like  that  of  experimental  shell 
number  23.  Fig.  14  was  plotted  from  the  analysis  of  "free"  deflection, 
that  is,  with  both  w^  and  b  as  independent  parameters.  The  calculations 
were  terminated  in  the  post-buckling  range  when  it  seemed  that  the 
pressure  was  no  longer  adequately  defined.  Fig.  15  was  for  the  same 
configuration,  but  for  a  case  of  volume  control.  Again  the  calculations 
were  terminated  at  a  point  in  the  large  deflection  range. 

A  number  of  analyses  were  made  for  the  "free"  deflection  case 
and  the  resulting  load-deflection  curves  are  plotted  in  Fig.  16. 

(The  buckling  pressures  for  these  analyses  were  listed  in  Appendix  d). 
From  this  figure  it  can  be  seen  that  the  buckling  loads  were  much 
too  high,  averaging  about  1.3  times  the  linear  buckling  load.  There 
is  some  variation  between  the  values  of  the  peak  pressures  and  there 
is  even  a  difference  for  the  two  cases  with  equal  X .  The  latter 
circumstance  points  up  the* fact  that  A  was  not  the  only  parameter 
needed  to  specify  shape;  o<  was  also  required  in  the  analysis.  This 
is  of  course  due  to  the  use  of  c  sin  cx  ^  sin  c  CX  and  c/b  sin  b  cX  , 
which  would  all  be  designated  as  cd  or  in  the  shallow  shell 
analysis.  Hence  it  is  quite  possible  that  X  is  not  the  only  impor¬ 
tant  parameter  for  deep  shell  buckling,  but  some  depth  measure  such 
as  (?<  or  H  is  needed. 
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Since  the  computed  buckling  loads  were  much  above  experimental 


results  and  even  above  the  calcr’ated  buckling  pressures  of  Budiensky, 
it  was  apparent  that  the  dexlected  shape  was  too  restricted.  Not  only 
was  the  only  form  of  deflection  a  central  dimple,  but  no  variations 
were  possible  for  the  mode  shapes. 

The  pressure-volume  change  curves  of  Fig.  14  and  15  show  that 
the  initial  portion  of  the  loading  curve  should  be  linear,  at  least 
for  large  X,  This  was  indeed  the  case  in  the  experimental  work 
(see  Fig,  lo). 

There  was  some  difficulty  in  finding  the  proper  range  of  b  for 
each  shell.  As  may  be  seen  in  Appendix  D,  the  values  of  b  for  the 
peak  of  the  curve  varied  in  no  orderly  fashion.  Further  it  was 
found  that  b  varied  rapidly  near  buckling  for  some  cases  and  very 
slowly  for  others.  At  least  the  values  of  b  which  were  of  interest 
were  always  near  90^  of  c.  In  all  cases  the  value  of  b  increased 
as  the  loading  proceeded,  reached  a  maximum  value  early  in  the  post- 
buckling  range  and  then  decreased.  Thus  most  of  the  sets  of  equations 
had  two  sets  of  solutions,  one  for  small  deflection  and  the  other  for 
large . 

The  difficulties  encountered  in  finding  equilibrium  positions  at 
large  deflections  have  already  been  discussed  in  the  Solution  Section. 
Evidently  more  nonlinear  terms  are  required  in  the  equation  for  77^ , 

This  would  mean  constructing  a  more  appropriate  order  of  magnitude 
criterion  than  use-j  herein.  Then  strain-displacement  relations  could 
be  simplified  using  the  new  criterion.  Probably  the  equations  of 
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Love  used  in  this  report  should  be  reviewed  to  determine  if  they  are 
valid  for  very  large  deflections. 

Some  of  the  added  problems  Just  described  may  be  avoided  or  sim¬ 
plified  if  the  loading  problem  is  divided  into  two  parts:  the  buckling 
problem  and  the  large  deflection  problem.  A;  least  for  the  buckling 
problem,  the  present  order  of  magnitude  criterion  and  strain-displace¬ 
ment  relations  are  quite  adequate. 

The  deflected  shape  used  in  this  analysis  is  not  recommended 
for  a  future  attempt.  A  reasonable  formula  for  deflection  should 
satisfy  all  boundary  conditions,  conform  to  the  known  and  expected 
mode  shapes,  have  as  few  teims  as  possible  per  parameter,  and  be 
reasonably  integrable. 
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CHAPTER  6  SUMMARY  AND  CONCLUSIONS 


In  the  experimental  portion  of  this  study,  thirty  nine  plastic 
spherical  shell  segments  were  tested  under  external  pressure.  Most 
of  the  tests  were  made  with  displacement  control  so  that  buckling 
occurred  without  a  change  in  average  displacement.  A  few  tests  were 
also  made  with  pressure  control  since  the  latter  situation  best 
approximates  a  practical  case.  In  this  group  of  shells,  the  shell 
shape  parameter,  A  ,  varied  from  3.5  to  25.  The  test  procedures  used 
in  the  experiments  were  listed  in  some  detail  because  of  the  important 
effect  of  creep  and  time  on  the  buckling  loads. 

All  the  shells  were  buckled  a  number  of  times  and  the  buckling 
loads  were  found  to  be  repeatable  under  unchanged  conditions. 

With  the  tests  under  displacement  control  a  measure  of  deflec¬ 
tion  was  known  at  each  stage  of  loading.  Hence  a  load  versus  deflec¬ 
tion  curve  was  plotted  in  some  cases.  Important  points  on  this  curve, 
were  the  buckling  pressure,  pressure  to  which  the  shell  Jumped  at 
buckling,  pressure  at  unbuckling  and  the  minimum  pressure  on  the 
curve.  The  critical  pressures  were  found  to  fall  in  a  fairly  narrow 
range  in  the  midst  of  test  results  of  other  investigators,  the  results 
for  constant  pressure  tests  being  slightly  above  those  with  volume 
control.  The  critical  pressures  exhibited  a  slight  downward  trend 
with  increase  in  A  and  the  minimum  pressure  was  evidently  a  function 
of  shell  depth  instead  of  A  •  The  other  two  important  pressures  were 
definitely  functions  of  A. 


91 


A  discussion  of  the  occurrence  of  buckling  is  presented  in  Section 
4D,  describing  the  variations  of  buckle  sizes  and  positions,  the 
shape  of  the  loading  curves,  and  the  changes  in  the  buckling  phenomenon 
with  changes  in  “A  . 

The  effect  of  imperfections  in  the  shell  were  investigated,  but 
no  relationship  between  imperfection  and  critical  load  was  found. 

The  accuracy  with  which  the  shell  was  clamped  in  the  test  chambers 
had  a  noticeable  effect  on  the  critical  load.  Since  the  plastic  of 
the  shells  exhibited  some  creep,  this  also  changed  the  critical 
loads.  A  simple  analysis  of  the  creep  effect  was  made  and  used  to 
cancel  out  this  effect. 

The  energy  criterion  of  Tsien  was  reviewed  and  compared  with 
experimental  results.  However,  Tsien  predicted  buckling  loads  under 
constant  pressure  would  be  about  half  those  under  constant  volume 
(displacement  control)  and  this  was  not  verified  by  the  experiments. 

The  tests  have  served  as  an  indication  of  the  value  of  plastic 
as  a  testing  material  and  of  the  importance  of  repeatable  testing  of 
the  same  specimen. 

A  variational  or  energy  method  was  used  in  the  analytical  study. 
The  equation  for  total  potential  energy  for  s3Tnmetric  buckling  in 
deep  spherical  shells  was  derived.  A  formula  for  deflected  shape 
was  chosun  based  on  the  shape  seen  in  the  experimental  work.  The 
equilibrium  equations  derived  from  *-he  energy  equation  were  found 
to  be  adequate  in  the  range  of  deflection  near  buckling  but  not  so 
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good  for  large  deflections.  For  an  analysis  far  into  the  post-buckling 
range  a  new  order  of  magnitude  analysis  and  new  strain-displacement 
relations  woulJ  probably  be  needed. 

The  equilibrium  equations  were  solved  to  find  equilibrium  positions 
and  thence  load-deflection  curves.  In  two  cases  a  plot  of  energy  versus 
deflection  was  also  made.  Analyses  were  made  for  five  shells  with 
from  14  to  25.  Since  these  were  considered  deep  shells,  it  was  found 
necessary  to  specify  both  A  and  oi  .  The  curves  were  quite  similar  to 
those  plotted  during  the  experimental  study  but  ctitical  loads  were 
3  or  4  times  too  high.  That  the  analytical  buckling  pressures  are  even 
above  those  of  other  Investigators  is  probably  due  to  a  restricted 
form  for  deflected  shape  which  failed  to  account  for  small  undulations 
which  occur  in  the  surface  prior  to  buckling. 

There  is  a  large  spread  between  theoretical  and  experimental 
buckling  pressures  and  even  quite  a  range  between  test  results.  This 
scatter  is  probably  due  to  clamping  conditions  during  testing  (5^ 
or  10^  but  possibly  a  much  larger  effect),  imperfections  in  the  shell 
(unknown  effect),  and  nonaxisymmetry  of  buckling.  The  tests  indicate 
that  the  latter  effect  may  be  rather  minor. 
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CHAPTER  7  SUGGESTIONS  FOR  FURTHER  STUDY 


This  investigation  has  developed  many  questions  and  problems 
which  still  need  a  solution  Among  these  are  the  following  topics 
for  study: 

1.  Unsymmetric  modes  of  deformation  should  be  considered  in 
determining  buckling  loads.  Techniques  with  a  reasonable  chance  of 
success  in  the  project  are  those  of  Weinitscke,  Thurston,  Budiansky, 
and  Caseman . 

2.  The  correlation  between  edge  restraint  and  buckling  load 
should  be  studied.  In  practical  situations,  the  edge  support  will 
not  be  clamped  or  pinned  but  merely  restrained  to  some  degree. 

3.  A  general  theory  for  the  effect  of  imperfections  on  shell 
buckling  is  much  needed.  First  the  mode  shapes  for  any  /\  must  be 
known  and  then  the  interaction  of  imperfection  size  and  position 
and  modal  deformations  could  be  considered. 

4.  The  equations  of  equilibrium  should  be  derived  and  a  study 
made  of  the  importance  of  terms  for  any  amount  of  deformation.  Many 
such  studies  have  been  made  but  usually  a  priori  and  not  based  on 
any  conrete  results, 

5.  It  would  be  useful  to  make  some  experiments  with  as  little 
error  caused  by  testing  conditions  as  possible  and  true  constant 
pressure  tests  should  be  performed  as  well  as  constant  volume  tests. 
To  reduce  the  effects  of  minor  imperfections,  the  shells  should  be 
thick  and  the  shells  should  be  clamped  between  rings  with  a  spherical 
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surface  so  that  no  rim  would  be  required. 


6.  The  effect  of  rate  of  change  of  load  over  the  surface  should 
be  considered.  So  far  uniform  loads  have  been  most  popular,  but  wind 
loads  and  point  load.s  also  produce  buckling. 
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FIO.  I  NOTATION  ON  SHELL  ELEMENT 
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FIG.  7  EFFECT  OF  SLOPE  ON  BUCKUNG  LOADS 
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FIG.  10  LOAD-VOLUME  CHANGE  CURVES 


i 


o 


o 

K> 


O  O 
CM  e 


•» 


0> 

N 


N 

m 


10 


112 


O 

U. 


BUCKLING  PRESSURE 
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FREE**  DEFLECTION  CURVES 


-w,+w,,  CENTER  OEFLECTION/raDIUS 


FIQ.  16  LOADING  CURVES  FROM  ANALYSIS 
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APPENDIX  A  EXPERIMENTAL  RESULTS 


The  data  for  Table  1  were  taken  from  tests  of  spherical  plastic 
shells  performed  ir  accordance  with  the  proceaures  described  in  the 
Experimental  Chapter  of  this  report.  Shells  to  be  tested  were  clamped 
into  the  testing  chamber  and  buckled  a  series  of  times.  Then  the  chamber 
was  opened  and  the  shell  was  reclamped.  With  the  slightly  altered 
boundary  conditioiis,  a  new  series  of  tests  was  run.  The  main  tests  were 
made  under  water  pressure  to  maintain  control  of  the  displacement.  The 
few  tests  made  under  air  pressure  were  so  designated.  When  the  same 
series  letter  was  applied  to  both  water  and  air  pressure  tests,  it  meant 
that  both  sets  of  tests  were  run  without  changing  the  clamping  conditions. 

The  buckling  loads  were  normalized  using  Eq.  (4-4)  to  make  the 

values  of  P^^  correspond  to  Ej .  Since  the  results  for  each  series  were 

nearly  identical  and  corresponded  to  one  set  of  conditions,  average 

values  of  (P  )„,  P, ,  and  P  were  recorded  for  each  series. 

'  cr'N  *  un 

The  buckle  and  imperfection  positions  are  given  as: 
c  center  of  shell 

n  near  center 

m  midway  between  center  and  edge 

e  next  to  edge  of  shell 

Following  the  values  of  q^,  X  .  a/h,  and  cX  for  each  shell  there  is  a 
brief  description  of  observed  surface  imperfections.  If  the  edge  strip 
or  rim  curled  up  (with  the  shell  convex  up)  after  forming,  a  positive 
slope  was  indicated.  An  estimate  of  the  height  of  bumps  and  their 
location  on  the  shell  was  made:  This  information  follows  the  slope. 
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In  the  table  E 


,  =  498,000  psl  and  g.  the  ic i sri't ,  is 


0.011. 
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Table  1  -  Buckling  Loads 


Series 

(p  ) 

'  cr^N 

P2 

P 

un 

Buckle  Position 

Shell 

#1 

=  .237 

psi  ,  "X  =  9 

.25,  a/h  = 

1635, 

=  7.40° 

L 

slope 

=  -.006,  imp.  =  .015", 

,  at  c 

) 

a 

.319 

.175 

.251 

n 

b 

.296 

.160 

.240 

n 

c 

.295 

.164 

.237 

n 

ave 

.303 

.167 

.243 

Shell 

#2 

^  .272 

psi  ,  X  =  9 

.20,  a/h  = 

1525, 

CX  =  7.62° 

1 

slope 

=  -.06,  imp 

.  =  none) 

a 

.265 

.256 

,256 

n 

b 

.253 

.243 

.246 

n 

c 

.240 

.226 

.229 

n 

ave 

.253 

.242 

.244 

Shell 

P 

(q^ 

=  .670 

psi,  X  =  12.3,  a/h  = 

972,  =  12.75° 

slope 

=  — .03 .  imp 

.  =  none ) 

a 

.400 

.167 

.261 

Cl 

b 

.390 

.156 

.264 

m 

ave 

.395 

.161 

.261 

d (air ) 

.368 

— 

— 

Shell 

#4 

(Qt 

=  ,78? 

psi  ,  X  ■=  4 

.67,  a/h  = 

898,  =  4.11° 

X 

slope  : 

=  /.02,  imp 

.  =  none ) 

a 

.373 

.363 

c 

b 

.359 

— 

— 

c 

c 

.348 

.342 

— 

c 

ave 

.360 

.352 

— 

S*-  11 

P 

(q^, 

=  .494 

psi ,  X  =  3 

,43,  a/h  = 

1135, 

=  3.29° 

slope  ; 

=  /.04 ,  imp 

.  =  none) 

a 

.610 

— 

b 

.698 

— 

— 

c 

c 

,556 

.545 

,545 

c 

ave 

.588 

,545 

.545 
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Shell  #7 

(q^  =2.275 

psi  ,  >.  =  4.92. 

a/h  =  527, 

cX  =  6. 

slope  = 

/  r\A  j _ 

f  .  ,  imp .  =  . 

ux  It  X  m  ; 

a 

.352 

— 

— 

m 

b 

.376 

.367 

— 

n 

c 

.353 

.345 

— 

e 

ave 

.360 

.356 

— 

Shell  #8 

(q.j.  =  2.375 

psi ,  X  =  5.15, 

a/h  =  516, 

cX  =  7. 

slope  = 

-.015,  imp.  = 

.0004"  at  m 

) 

a 

.313 

.307 

— 

n 

b 

.315 

— 

— 

n 

c 

.294 

.290 

— 

n 

ave 

.307 

.298 

— 

Shell  #9 

(q^  =  1.320 

psi ,  =  7.32, 

a/h  =  692, 

=  9. 

slope  = 

none,  imp.  at  : 

m) 

a 

.447 

.216 

.375 

n 

b 

.475 

.215 

.380 

n 

c 

.491 

.199 

.371 

n 

ave 

.471 

.210 

.375 

Shell  #10 

(q^  =  1,400, 

P^=  7.32,  a/h 

=  672,  cX  = 

9.19° 

slope  = 

-.003,  imp.  none) 

a 

.449 

.229 

.381 

n 

b 

.442 

.217 

,372 

m 

ave 

.445 

.223 

.376 

Shell  #11 

{q^  =  1.435, 

A  =  7.25,  a/h 

=  662,  cX  = 

9.10° 

slope  = 

-.003,  imp.  = 

.001”  in  symmetric 

a 

.447 

.205 

.366 

m 

b 

.436 

.201 

.358 

m 

ave 

.442 

.203 

.362 

Shell  #12 

{q_  =  1.350, 

X=  7.47,  a/h 

=  684 ,  cx  = 

9.24^ 

no  slope,  imp.  none) 

a 

.452 

.211 

,355 

— 

b 

.466 

.211 

,371 

ave 

.459 

.211 

.363 

Shell  #13 

{q^  =  1.352, 

X  =  7.34,  a/h 

=  693,  cX  = 

9.08® 

no  slope 

,  central  hollow) 

a 

.456 

.213 

.388 

b 

.459 

.210 

.381 

ave 

.458 

.212 

.385 
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Shell 


Shell 


#14 

(q* 

a  .189  psi,  %=  5.2, 

a/h  =  1735,  o(  =  4.03° 

X 

slope  =  /.02,  Imp.  = 

.01"  at  m,  1"  dia.  bump) 

a 

.428  .333 

.404  c 

c 

.485  .315 

m 

c 

.472  .338 

.432  m 

ave 

.462  .329 

.418 

#15 

(q* 

=  .209  psi,  X=  5.1, 

a/h  =  1740,  eV  =  3.94° 

1 

slope  =  /.02,  point 

imp. ) 

a 

.544  .377 

.467  n 

b 

.515  .396 

.462  m 

c 

.494  .323 

.428  m 

ave 

.518  .365 

.452 

Shell 


#16 

(q^  =  -663  psi 

,  7.03,  a/h  = 

977  ,  0<  =  7 . 27° 

no  slope. 

imp.  =  .023"  at  m 

in  symmetric  ring) 

a 

.425 

.269  .362 

e 

b 

.468 

.263  .381 

n 

c 

.460 

.260  .381 

n 

ave 

.451 

.264  .375 

Shell 


Shell 


Shell 


Shell 


#17 

(q^ 

=  .629  psi,  =  6.87, 

a/h  = 

1000,  o<  =  7.02° 

no  slope,  point  imperfection  at  c) 

a 

.493  .238 

.388 

n 

b 

.496  .240 

.372 

m 

c 

.480  .240 

.387 

m 

ave 

.490  .239 

.382 

#18 

(q- 

a  .638  psi,  >.  =  6.86, 

a/h  = 

995,  Ok  =  7.02° 

1 

no  slope,  imp.  =  .01" 

at  e. 

various  bumps) 

a 

.535  .253 

.417 

n 

b 

.500  .238 

.402 

n 

c 

.493  .229 

.389 

n 

ave 

.509  .240 

.403 

#19 

(q™ 

=  .518  psi,  5.76, 

a/h  = 

1105,  o<  =  6.57° 

1 

slope  =  /.03,  hollow 

at  c) 

a 

.447  .272 

.394 

m 

b 

.515  .269 

.405 

m 

c 

.490  .292 

.417 

m 

ave 

.484  .278 

.405 

#20 

(q-. 

=  .522  psi,  X-  3*31. 

a/h  = 

1100,  CX  =  3.32° 

1 

slope  =  /0.4,  imp.  = 

.01"  at  n) 

a 

.549  .549 

— 

c 

b 

.571  .571 

— 

c 

ave 

.560  .560 
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Shell 

#21 

(q_  =  .582  psi 

,  X  =  12.0,  a/h  =  1040, 

— 

r 

1  O 

no  slope. 

imp.  =  .013"  at  m) 

a 

.514 

.155  .302 

n 

b 

.512 

.154  .296 

m 

c 

.533 

.187  .296 

- 

ave 

.520 

.165  ,298 

Shell 

#22 

(q„  =  1.620  psi,  X=  16.1,  a/h  =  625, 

ok  = 

21.3° 

slope  =  -. 

015,  imp.  ,01"  at  c) 

a 

.325 

.148  .230 

n 

b 

.317 

.147  .226 

n 

c 

.330 

.153  .225 

n 

ave 

.324 

.149  .227 

Shell 

#23 

(q  =  1.535  psi,  X=  16.1,  a/h  =  642, 

cx  = 

20.9° 

no  slope. 

imp.  =  .01"  at  c) 

a 

.459 

. 135  . 236 

n 

b 

.463 

.143  .235 

n 

c 

.460 

.142  .213 

n 

ave 

.461 

.140  .238 

Shell 

#24 

(q„  =  .495  psi 

,  X  =  9.85,  a/h  =  1130,  C<  = 

9.48° 

slope 

=-.015,  imp. =.02"  at  e) 

a 

.525 

.184  .270 

c 

b 

.472 

.187  .338 

m 

c 

.492 

.187  ,350 

n 

ave 

,496 

.186  .319 

Shell 

#25 

{q,j,  =  .481  psi 

,  X  =  9.9o,  a/h  =  1145, 

cx  = 

9.48° 

no  slope. 

hollow  at  c  of  unknown  magnitude) 

a 

.259 

. 234  . 238 

c 

b 

.313 

.271  .296 

c 

ave 

.286 

.252  .267 

Shell 

#26 

(q^  =  .221  psi 

,  A  =  11.6,  a/h  =  1695, 

cx  = 

9.10° 

no  slope. 

imp.  =  .006"  at  m) 

a 

.328 

.187  .271 

c 

b 

.291 

.157  .174 

c 

ave 

.309 

.172  .222 

Shell 

#27 

(q,j.  =  .228  psi 

,  A  =  11.7,  a/h  =  1670, 

cX  = 

9.24° 

no  slope. 

imp.  =  .003"  at  e) 

a 

.421 

.228  .340 

e 

b 

.446 

.205  ,350 

e 

c 

.404 

.169  .314 

e 

ave 

,424 

.201  .335 
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Shell  #33  (q^ 

=  .616 

psi  ,  X  =  20.2  , 

a/h  = 

1015, 

CX  =  20.9° 

no  slope,  imp.  =  .015 

at  c 

,  rise 

) 

a 

.339 

.115 

.192 

b 

.450 

.104 

.203 

n 

ave 

.395 

.109 

.197 

b(air) 

.432 

— 

— 

e 

Shell  #34  iq 

-  .606 

psi  ,  X  =  20.3 , 

a/h  = 

1020, 

CX  =  20.9° 

uo  slope,  imp.  =  .02" 

at  c , 

rise) 

a 

.319 

.109 

.184 

— 

b 

.421 

.099 

.207 

n 

c 

.409 

.100 

.182 

n 

ave 

.383 

.103 

.191 

b(air) 

.415 

— 

— 

n 

c (air ) 

.411 

— 

— 

n 

Shell  #35  (q^ 

=  .243 

psi ,  X  =  25.6, 

a/h  = 

1615, o'  =  20.9° 

no  slope,  imp.  =  .02" 

at  c , 

rise) 

a 

.277 

.033 

.116 

n 

b 

.283 

— 

e 

b 

.329 

.107 

.181 

n 

c 

.331 

.105 

.182 

n 

ave 

.305 

.082 

.160 

Shell  #36  (q^ 

=  .243 

psi  ,  X  =  25.6, 

a/h  = 

1615,  CK  =  20.9° 

no  slope,  imp.  =  .02" 

at  c. 

rise) 

a 

.259 

.054 

.135 

e 

b 

.310 

.042 

.141 

e 

c 

.293 

.037 

.124 

n 

ave 

.287 

.045 

.133 

a(air) 

.315 

— 

e 

b(air) 

.349 

— 

— 

e 

Shell  #37  (q^ 

=  .272 

psi,X=  24.7,  a 

/h  =  1525,  W  =  20.8° 

no  slope,  imp.  =  .01" 

at  c , 

rise ) 

a 

.406 

.119 

.206 

e 

b 

.335 

.121 

.171 

e 

c 

.352 

.082 

.139 

e 

ave 

.373 

.107 

.172 

a(air) 

.390 

— 

— 

e 

^38 

(q,j.  =  .600 

psi  .  X=  12.4,  a/h  = 

1025,  0(  =  12.6 

slope  : 

=  -oil,  imp.  at  e) 

a 

.510 

.148  .262 

m 

b 

.589 

.150  .274 

m 

c 

.557 

.  165  . 300 

— 

ave 

.552 

.154  .279 
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V 

Shell 

#39 

(q™. 

=  .700  psl,  X  =  12.6, 

a/h 

950,  O'.  =  13,3° 

V 

1 

slope  =  -.006,  no  imp, 

.) 

a 

.309  .116 

,168 

m 

b 

.267  .128 

.170 

- 

ave 

.288  .122 

.169 

Shell 

#40 

{q^ 

=  .616  psi ,  X  =  12.4, 

a/h  = 

1010,  C5i=  12.6° 

< 

no  slope,  no  imp.) 

a 

.425  .144 

.209 

e 

b 

.436  .137 

.205 

m 

c 

.440  .130 

.197 

~ 

ave 

.434  .137 

.204 
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APPEK'IUX  B 


Modulus  and  Creep  Coefficient  Determination 

Beam  Teats 

Samples  of  the  shell  material,  cut  to  6‘*  by  1^**,  were  tested 
as  simple  beams  to  determine  the  modulus  of  elasticity.  The  sample 
numbeirs  correspond  to  the  shell  nximbers  since  they  came  from  the 
same  sheet. 


Table  2  - 

Modulus  frcan 

Beam 

Tests 

No. 

Thickness  Width 

psi 

No. 

Thickness 

Width 

El,  psi 

U 

.oaae** 

1.UU5" 

1*95,000 

12 

.0501*" 

1.1*90" 

503,000 

5 

.0830 

l.UUd 

500,000 

13 

.0500 

1.1*56 

500,000 

6 

.0777 

1.U56 

1*89,000 

lU 

.oia9 

1.515 

509,000 

7 

.0829 

1.U29 

1*77,000 

15 

.01*27 

1.1*80 

1*91*,000 

8 

*0783 

I.U46 

532,000 

16 

.0^29 

1.519 

1*90,000 

9 

.0512 

1.U12 

1*91,000 

17 

.0la9 

1*515 

1*98,000 

10 

.0508 

1.510 

1*92,000 

18 

•01*21* 

1.1*37 

502,000 

11 

.0508 

1.1*77 

1*79,000 

19 

.oia9 

1.1*20 

521*,C00 

The  average  value  of  is  U98,000  psi  and  the  standard  deviation 
of  this  value  ia  l.W 
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Tensile  Test;,; 

The  tensile  tests  were  nin  using  SR-U  strain  gages  and  a 
dead  weight.  The  gages  and  the  glue  attaching  them  to  the  specimens 
take  a  certain  amount  of  load.  To  take  account  of  the  added 
stiffhess,  the  apparent  modulus  can  be  modified  according  to  equation 

®1  “  ®apparent^^  “  area  ^ 

in  which  B  describes  the  gage  stiffhess.  If  B  and  are  assumed 
to  be  conatants  for  the  teneile  specimens,  a  least  squares  formula¬ 
tion  can  be  used  to  give 

-  5lU,000  psi 

B  -  .00687 

The  standard  deviation  of  the  value  for  is  1.8j{. 


Table  3 

-  Modulus 

from  Tensile 

Tests 

Ho. 

Thickness 

Width 

E 

apparent 

®1  “  ®app^^  ^ 

5 

.0829 

1.UU7 

556,000 

52U,000  psi 

6 

.0789 

1.U57 

5U5,000 

512,000 

8 

.0770 

Uhh9 

539,000 

506,000 

9 

.0512 

i.iau 

578,000 

523,000 

16 

.01429 

1.528 

556,000 

U98,000 

17 

•OU26 

I.I462 

583,000 

518,000 

wt 


3  L'S 


The  adjustment  for  rage  effect  is  from  6  to  11/J  and  suggests 
that  this  tensile  modulus  is  rather  questionable. 

Gages  placed  crosswise  on  some  of  the  specimens  gave  a  value 
for  Poisson's  ratio  of  0«iil* 

Creep  Coefficient 

Ihe  curve  of  strain  versus  time  for  the  tensile  specimens 
follows  the  equation 

fairly  well  in  tt»  range  of  interest.  Denoting  the  strain  at  one 
minute  by  ,  and  that  at  10  minutes  by  ^/o  ,  the  value  of  g 
can  be  foxmd  from 


The  tensile  tests  then  give  the  values  found  in  Table  k 


Table  ii  -  Creep  Coefficient 


No. 

€/p/e, 

E 

Mo. 

g 

5 

1*022 

•00955 

9 

1.021 

.00911 

6 

1.020 

.00667 

16 

1.039 

.01690 

8 

1.029 

.01255 

17 

1.026 

.01128 

The 

average  value  of  g  is 

.0113. 
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B-'j-rm  End  Rotation 


Since  the  simply  supported  beam  used  in  the  moduiue  deter¬ 
mination  undergoes  considerable  deflection,  the  end  rotation  jiwiy  ba 
a  nonlinear  function  of  load.  To  determine  this  relationship,  let 
P  *  load  on  the  beam 

i  -  length  of  beam  he^t^/eert  iijtec/ 

M  ■  moment  in  beam 

^  "  rotation,  rotation  at  ends,  rotation  in  the 
linear  theory 

X  "  variable  horizontal  distance  measured  from  center 
a  >  variable  distance  along  beam  from  center 
8^  “  value  of  s  at  X  ^ 
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for  the 


In  the  nonlinear  term  let  ^  ** 

a 

integration  and  ^  evalur.tion.  Using  these 

approximations  in  the  equation  for  x  gives 


/ 


*  -y 
£  ^ 


Sa  - 


2VO 


Therefore 


i/« 


but 


y  - 


Then  the  first  nonlineer  approxlnetion  is 


n. 
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APPENDIX  C  ENERGY  EQUATION 


The  energy  equation  is 


IT  =  E'o^^ 


(5-34) 


This  equation  is  obtained  by  performing  the  integrations  of  Eq.  (5-25), 
rearranging,  and  eliminating  terms  which  cancel  with  the  aid  of 
Eq.  (5-29). 

If  the  integrals  of  Eq.  (5-26)  were  evaluated  exactly,  many 
terms  would  be  obtained  which  are  as  small  as  those  neglected  already. 
Retaining  such  small  terms  is  unreasonable  and  so  a  criterion  was 
set  up  to  eliminate  th«m..  The  evaluated  energy  function  contains  the 


following  types  of  terms: 


-T-yj  arTd  - €■ 

b 

These  two  types  are  basically  different  since  the  exponeiktlal  is  not 


readily  compared  with  l/b  and  also  because  the  exponential  term  increases 
on  differentiation  while  the  other  decreases.  Therefore  the  highest 
order  terms  of  each  type  were  retained.  This  criterion  allowed  for 
liberal  use  of  approximations  in  evaluating  since  this  whole 
term  is  essentially  of  secoisd  order. 

The  E's  of  the  equation  will  be  defined  below  in  the  same 


manner  as  they  were  for  the  Fortran  statement  used  in  their  evaluation. 
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Bj,  =  /  6^  "  ^  8 

B. ^/.9s'?  C^  -J?.^ir?  sj,to<ii7^oi  7*  /i 

Bof/,^:r7  C^  '/.479r/r,M  Ar^ct  >5 

C^  7- 

j  /  ^dr'7  s/pioi  f%  /5a 

I B^7  - /.-^s^ s‘/,fttn  ^3 

The  G’s  and  P's  are  given  by  the  following: 

C,  =S/?,(K3,-^lf„)-i.‘ff/,(l^:,*^Kk)  -f 

¥  -COSot)-^  S^ 

^  ^ - g^X ^  yA^ C^9M  ^3^ 

^  B/ •“  COS'O^  ¥"  -^ 

^  **  C^^j  ¥  ^3  )  ¥  ^,  ^  B/ ^  C^S*»t)  “■  yM^3  ’“CcS^J 

=(^33  '  ^3^X f Ccrc^)  ¥  S^ 


r.  = 


(  / 


Cs//fOi  J?c\y}r^<^ 


>5  ^  ^.sjfcosBai  fSf'rtlM  ¥^CasA>^ca^o/J&'^  -  -  _  /i^C^So^  ^ 

^x//r  U  Jt 
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p  .^3? 

*3  "  c  ““  ch^  '  'IP' 


2.  ^:i^o( 


^  Si»  k>e{  -^S^SC-^A)  C0xlot  ^  c  ~:>6C0X^4*>t  y  ^M<r  Dj^o/1 

y, asa 7S‘/S^  )y  J}kA  y  \  „  2^3  /  ^g^jco/o/  l^^CaSM 

»  c /*  CoS»i  {/^COFaH.)  J  C*  /yrojro^ 

^{.v/?^ ^  cosot  ^ j^i  ff.^j:p£^.  ) 

^J>./^9 A^  A^/-^  y_JL  y  ^  cosa  -  ,  /^S'cor^^i  /  rorw 

d  ^  <3»rw  /yoflr«« 

-  yycasoi  )  -  3  S/O  *oi  -^,9 CoS  ex  -./COXC)} 

St  /  Cos  ot 

**  CCS  Ot) 


f^  =  sAtoi  ^AcosSAc(S//lct  ¥As/cJtAmt 


ry/fci 


^COSeat  ^-^COS^eoTot  ¥  "9 COSat  ■‘^S/zpAm.  ’-^cosAei 
'•S/oSUd  -  cvsSiAot)  -  “  9j>/r^  (f“^*^sy\fAd  s/m  Sxt 

y*  y  cosA*i)cos  V  y “"  cosot) 


^  5=  ,  0?S'c  Of/oi  y  /  ^3/ 

/  y  coso>( 


-.?  //^ccsoi{y  -  ^ 


/  y  cosot 


caSfi 
S/ri^ 


^  V 


^ -/  -  -g^r!  y  ccssAoi  oi^y 


1 


CO^  2 ^oi 


ccx^oi ^ ^  i^.^COS3^  vx/x/ir 

ya-  yF' 


^ /^p  ^  sy/td  —•  9^cas‘^o»  tyy>f^*>i)c0Stij  -^±7^3  cosei, 

6^syh\( 


^J!6»^/o.?-:r.ysyyr  It  »>  £.^±'?S££^  )  i^ 

/^  ‘  /><  CoS  •I  / 


/  i^Cosoi 


Ihe  H's  and  R's  and  S' a  are  in  turn  defined  as  follovst 


=:  ,  COSOK 

C  -Jc*" 

ev>»  cK 


// - L _ )^ 

'•  Ses/j0Uy 


cosu 


*  C^rAl 


=* 

=  Csypt^u  -  /t^CcSo^ 


^yyiV  -  xvyfr  Xrf  <J>7^  y  cosJ^  cc 

7~~ - :  /-/  /  coriot  -  r/s^co^tJ  -y  cosMcd/a^ 

^^C/’-caro/)L 


y' .  (^/'/^CasoO 


R..  = 


C  s/^  U  -  COS  CK 
<r4i>so^J 


- 


cA  S^/fOi  -{'Cf^dJ s//7l^ 

'"-COSO*  I 

^ -  s  ^  ^ ^  cofc/  —  "y 

<rO-W0f)J 

^  ^jfO^Cafofj^X/W^  ccrlmf  J 


^^Aisl. 

/-CaS'tA 


^3/ 


^O'-  CaSai) 

t* _ ^  / •§■  -CoSoi  ^  cosit  ) 

S^/  -  trasof)  3  X 


*/9X^ coso^  /j  _  ^Cc/ot ) 

/—  C^aSot  '  ' 


Tjjf  a 


.?cA  ai>Soiff~^f-^i cosJot  r*  ^C^-AVs//r^^/CA/^s/lliJ7 
'-C.SO,  /  *  y 


/*  -  ^  r<pxo< 

/-  cosot 


^  /^  COS^Aot  1*-  ^AA  COsAmt  »* 

/-  COS  ot  « 7sr 

^^^sy*^Axco/<x  xJ^ccS-A^ccA*4  ’‘^r/>rjAoiC^  ■^^COS'^C^^l) 

y-~CoSe>t  6Y/-COS*Jt)  O^O-^COSUJ 

“XZ/jAet/COSJ^  ^  CosAft  CCS 3ol)  -  /coS^  -coslic/ 

^  /-CdX«i  t  f/ 

-  <afcrf  cr  ^ ^  A'cA  Sj/rAoi  -  ^oov^o' )  ^ 


Jy  = 
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-/  /Zy/  c  ^  y'>fAo^  7^  (s^ :ZZc.psAki  :ts.  t/j-^'-A 

^C  '*■-  -?«r^  * 

tZ2^3)o>sbu  -(3c  ^PVc/A^t  ^  "^^ZcyA^ccLtA^i  -  ^c  -  ^  Jj  ✓//  iZl 
■^/.^S'AZ  Cos  oi  -  .  //^  coZoi  Z^syA  Jiift  —  ) 

A 'f  CoSoi  /^/  -  Aco/o(  ) 

S*c. 

-  -/.^CA^c* Z?c-9AJcosJ^  C^c-^AJ r/A^ti 

L  /Xzic^-  Sc  A  ^  a^A 

^  ^  yc-9A)  s^J/7Ao(  y.  Z C-^lJcos-^A^  ^  Zc^JUJfyADAj 

f*  Z'^os^Ao/  r//fCfAofj] yc  ^  y9/Zf  e  Sei/Zc-JJcorAei'^  Zci^AjsArJai 

A  ^///,  Zcos-Aoi  A.dA/sS' /Zy  Y'^ZfZ  -  fA  'Sa.  y.  S'Cof  Sojj 

yfcAcos^&''^‘^f  Cc^yScA  ^:3AyaKrJ>A^^Zc^^ScA^aA^)sAt^U 

Z c  ^yA  VA  \) 

y-  /Zr  /Z,/AZji  Zees' -_  /Z// 


c\>A 


y  Z 9A  Aj  Xy^y/A^of  <?"  '^°^fZc^'-cA  ^A  \)cosAwi  y  c  '*^sA*9AJ7 
c^^A*'  ^  -Z 


'JX^IX  -  .  y^SfS"  <r  /^  /  COSJtAef  tS>  '*taIoi  f  cc/oi 

A?  ‘ 

^CoZa^  -  -  e'^Z^As-ye^fd!  CosZoi  - 


V* — ~  cos 
J?A 


^  -  CeysJoi )  V  (y  '^^}.C.o£Ai!i  T*  c  y/'i^U 


t{^c -PAJccs^^ -{Sc *?i)rj/fJ>A  ^  Cc^^jcctlj 

Coxot  -  ./J?2ir //jf  C€f/t3t :>t»i  - 
-/  //g^  CcSU  ^/  -  Sctf/oi  J 

^  3c A ar^  V 

/  ^^yc-9A)fj/tloi  ^  ^ ^  ^^JfygaL/ 

f*  ^  y9/^e  Jtj/tc-^JJcapjmt^  fc^d)s^Md 

/J?  y  ^ 

i'co^^oi  y  y-. y  S<or  Jc^ 

V,  y9c^c»Sc>i6!~^‘^fCc^^3c^  '^3^^€oij-J>4«e - y^c*‘-Sc^ 

Ccor^M rJS£iJ/!u 


^//yy^^of  e~ VcosAwi  y  «•  V>>»^<7 


—  cZ^^T”  -  .  /:fJfS'C’'^  r'rfetf  f  ^  car^^i  ^jnr)fJtJc4  y  :<Wa( 


y*/ ^  ^  :>cosS€^ 

9. 9f  9J>  cosl^  ^S/,f/>f^)x/>f'ii  ^{.3i/>f^>i  ' 

^S  3ot^  COS^  XZ/f^  /  ^{cz>S  '  S'y n  el  -  Coi3ol^l 

9.030i>JiS‘ CCSjL(  X^ec-teA 

■?“ 

■f£cos3^  -./a.*-^  ,*  ^)  ■/■  /.9i  Mj,e'^‘'7’-f^eos3^>i-/kM 

t  COS' ^eZl  M  jf^it  3 *>i  f§  COf-^eM  Cot 


-liiol 


'390 


4  (si  C0timiSM3ol^  ^  kJ  COT  Sol  -  isz'f}itl  s/if  SeJ  ^-pCOXiA 


XSicOsUdSMScl  ^  Si y/rfSSoiJ/H  aw  V  ifcOSSimt 


In  all  of  the  ebove^  Poiiieon's  ratio  hae  been  tak«)  aa  0*3» 
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APPENDIX  D  ANALYTICAL  RESULTS 


Th«  analyses  outlined  In  the  Solution  Section  waa  performed  on 
five  shell  configurations.  The  results  of  Table  S  are  for  buckling 
loads  which  were  obtained.  The  deflection  anl  the  value  of  b  at 
buckling  were  also  listed. 


Table  5 


P 

cr 


18. i 

20.9® 

1.368 

16.1 

15.0® 

1.27 

20.0 

26.0® 

1.372 

25.0 

20.9® 

1.31 

14.0 

18.0® 

1.20 

*1 

b/c 

000  37 

.001  31 

.83 

000  273 

.000  73 

.87 

000  413 

.001  35 

.97 

000  085 

.000  46 

.92 

000  45 

.001  27 

.93 
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